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ABSTRACT 

The  available  viscosity  and  thermal  conductivity  data  in  helium, 
argon,  hydrogen,  nitrogen,  are  reviewed.  Deviation  charts  vs  tem¬ 
perature  are  drawn  up  for  each  of  the  gases  of  interest  in  the  tem- 
per^^ture  range  200-2500  K.  High  temperature  thermal  conductivity 
irie.ajS’urements  in  air,  a!rgon  and  nitrogen  are  also  reviewed  at  tem¬ 
peratures  between  2000  and  12,  000  K.  The  mnethods  of  measuring 
thermal  conductivity  and  viscosity  are  discussed  with  emphasis  on 
the  extension  of  available  techniques  to  higher  temperatures  and/or 
pressures.  Major  emphasis  is  placed  on  constricted  arc,  heat  trans¬ 
fer,  schlieren,  and  ultrasonic  techniques.  Finally,  measurements  of 
helium  thermal  conductivity  in  the  temperature  range  ~400-3000  K  by 
the  ultrasonic  technique  in  the  shock  tube  are  presented. 


.  -.-jyn'': 


iii 


AEDC-TR.68.105 


CONTENTS 

FOREWORD  ii 

ABSTRACT  iii 

NOMENCLATURE  viii 

l.  INTRODUCTION  1 

n.  DATA  REVIEW  5 

2.  I  Deviation  Charts  5 

2.  2  Helium  6 

2.  3  Argon  8 

2.  4  Polyatomic  Gases  9 

2.5  Hydrogen  11 

2.  6  Nitrogen  12 

2.7  Air  14 

m.  METHODS  17 

3.  1  Classical  Techniques  17 

3.  2  Cascade  (Maecker)  Arc  23 

3.  3  Heat  Transfer  29 

3.  4  Schlieren  Technique  34 

3.  5  Ultrasonic  Technique  36 

3.  6  Ultrasonic  Measurements  of  Thermodynamic 

Properties  41 

3.  7  Laser  Acoustics  42 

3.  8  Ultrasonic  Absorption  Measurements 

in  Helium  47 

REFERENCES  50 

APPENDIXES 

I.  2^4ATHEMATICAL  RELATIONS  FOR  THE  END 

WALL  BOUNDARY  LAYER  63 


V 


AEDC-TR.68.105 


II.  ILLUSTRATIONS 

Figure  Page 

1  Viscosity  deviation  chart  for  helium  70’ 

2  Thermal  conductivity  deviation  chart  for  helium  71 

3  Excess  viscosity  for  helium  72 

4  Excess  thermal  conductivity  for  helium  73 

5  Viscosity  deviation  chart  for  argon  74 

6  Thermal  conductivity  deviation  chart  for  argon  75 

7  Thermal  conductivity  for  argon  in  the  range 

3000°K-12,  000°K  76 

8  Excess  viscosity  for  argon  77 

9  Excess  thermal  conductivity  for  argon  78 

10  Eucken  factor  for  monatomic  gases  as  a  function  of 

temperature  79 

11  Viscosity  deviation  chart  for  hydrogen  80 

12  Thermal  conductivity  deviation  chart  for  hydrogen  81 

13  Eucken  factor  for  hydrogen  82 

14  Viscosity  deviation  chart  for  nitrogen  83 

15  Excess  viscosity  for  nitrogen  84 

16  Thermal  conductivity  deviation  chart  for  nitrogen  85 

17  Eucken  factor  for  nitrogen  86 

18  High  temperature  thermal  conductivity  deviation  chart 

for  nitrogen  87 

19  Viscosity  deviation  chart  for  air  88 

20  Excess  viscosity  for  air  89 


AEDC.TR.68.105 


Figure 

I 

21 

T 

22 

\  '■ 

23 

’•  I 

24 

25 

r* 

I 

26 
27 

28. 

✓  * 

29 

30' 

O') 

31 

32 

*  r  j 

33 

F  I 

34 
l-H 

35 

•  ,) 

’  •’ > 

36 

•  / 

I  • 

37 

sa" 


Page 


Thermal  conductivity  deviation  chart  for  air  90 

Excess  thermal  conductivity  for  air  91 

High  temperature  thermal  conductivity  deviation  chart  for  air  92 

Cutaway  schematic  of  a  high  temperature  capillary  viscometer  93 

Thin  wire  thermal  conductivity  cell  94 

The  constricted  arc  95 

The  effect  of  arc  current  and  diameter  of  the  effective 

thermal  conductivity  96 

Heat  transfer  in  the  reflected  shock  region  97 

Experimental  apparatus  for  schlieren  measurements  98 

Schematic  of  schlieren  system  99 

Schlieren  photograph  of  end  wall  boundary  layer  100 

Test  chamber  for  high  pressure  sound  speed  measurements  101 

Velocity  of  sound  in  nitrogen  gas  at  room  temperature  102 

Velocity  of  sound  in  nitrogen  gas  as  a  function  of  pressure  103 

Scattered  intensity  spectrum  104 

Model  laser  gas  scattering  system  105 

Diagram  of  apparatus  for  ultrasonic  measurements  in  the 

shock  tube  106 

Experimental  viscosities  for  helium  107 


AEDC-TR-6B-105 


NOMENCLATURE 


Experimental  viscosity  and  thermal  conductivity 

T7  ,  X  Reference  viscosity  from  some  readily  available  tabula- 

IT  r 

tion.  Reference  values  are  tabulated  on  deviation  charts 


P 

f 

D, 

int 

T 

c 

T. 

1 

Z. 

1 

R 

P 

T 

D.. 

ij 


Percent  deviation  between  theory  and  experiment  (Eq.  1) 

Low  density  (~1  amagat)  viscosity  and  thermal  conduc¬ 
tivity 

Excess  viscosity  and  thermal  conductivity  over  the  low 
density  values 

Density 

Eucken  factor  defined  by  Eq.  (5) 

The  diffusion  coefficient  for  internal  energy.  In  some 
cases  internal  energy  may  be  passed  from  molecule  to 
molecule  faster  than  the  molecules  themselves  diffuse 

Mean  time  between  collisions,  taken  as  T  =  -7 

c  4  p 

Relaxation  time  for  ihe  ith  internal  mode  defined  by 
Eq,  (48) 

T. 

Collision  number  for  the  ith  mode  Z.  =  — 

1  T 

c 

Perfect  gas  constant  usually  used  in  dimensionless 
equations 

Pressure 

Temperature 

Diffusion  coefficient  of  the  ith  species  due  to  a  concen¬ 
tration  gradient  in  the  jth  species.  For  i  j  D^j  is 
the  binary  diffusion  coefficient  for  i  =  j,  the  self¬ 
diffusion  coefficient 
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Mole  fraction  of  the  i  species 

Capillary  radius  at  any  point  along  its  length 

Capillary  length 

Rayleigh  nvimber 

Grassoff  number  (see  p  20) 
rj  C 

Prandtl  number  - R- 


Heat  flux 
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Radiative  heat  flux 

i 

Total  radiative  intensity  from  a  cord  of  an 
_ erg  _ 


arc 


sec  cm  d£2 


Solid  angle  subtended  by  the  receiving  volume 
at  the  center  of  the  radiating  volume  element 


I 
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B 
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Radiative  intensity 


Planck  function 


a 
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•  Spectral  absorption  coefficient  for  light 


V 


Light  frequency 


T 

V 


Optical  path 


AT 

w 


End  wall  temperature  jump  behind  reflected 
shock  boundary  layer 


n 


Similarity  parameter  which  replaces  distance 
in  the  compressible  case  (Eq.  34) 
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Enthalpy 

Enthalpy  in  reflected  shock  region 

Enthalpy  at  the  wall  in  an  end  wall  reflected 
shock  boundary  layer 

Contrast  in  a  photographic  negative 

Specific  heat  (cps)  at  constant  pressure 

Specific  heat  at  constant  volume 

Contribution  of  the  ith  interval  mode  to  the 
specific  heat 


Specific  heat  ratio 

V 

Sound  frequency 
Thermal  diffusion  ratio 
Mole  fraction  of  the  ith  species 
Atomic  mass  of  the  ith  species 
Sound  speed  at  frequency  w 
Sound  speed  at  co  =  0 


Compressibility  factor  z  = 


pRT 


Wave  number  of  a  phonon  in  the  fluctuation 
spectrum 

Wave  number  of  incident  and  scattered 
light,  respectively 

The  instantaneous  fluctuation  in  density  at 
the  point  "?  and  the  time  t 
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Ztt 
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6(r) 
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Frequency  shift  of  scattered  light 

Boltzmann' s  constant 

The  delta  function 

Effective  collision  frequency 

Half  width  of  Brillouin  line  and  Rayleigh  lines, 
respectively 


y  Parameter  describing  the  ratio  of  phonon 

wavelength  to  mean  free  path  of  the  gas 
molecules 
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SECTION  I 
INTRODUCTION 


The  objective  of  this  work  is  the  experimental  determination  of 
gas  properties  at  high  temperatures  and/or  pressures,  with  emphasis 
on  transport  properties,  ^e  immediate  goals  of  the  present  work  are 
to  extend  existing  tabulations  for  immediate  engineering  needs  and  to 
provide  long  range  support  for  theoretical  work.  The  temperature 
range  is  300°K  to  20,  000°K,  and  the  pressure  range  is  1  to  10^  atmo¬ 
spheres.  The  principal  transport  properties  of  interest  are  viscosity, 
thermal  conductivity  and  diffusion.  In  addition,  consideration  of  other 
properties  such  as  electrical  conductivity,  optical  depth  and  certain 
thermodynamic  properties  are  of  interest. 

The  theory  of  energy  transport  in  dilute  gases  is  well  understood 
at  s\ifficiently  low  temperatures  that  chemical  reactions  may  be  neg¬ 
lected,  and  at  sufficiently  low  pressures  that  only  binary  collisions 
must  be  considered  (Refs.  1,2,3).  Two  types  of  practical  information 
are  provided  by  kinetic  theory.  First,  given  the  interaction  potential 
for  the  chemical  species  in  the  gas  of  interest  (Ref.  4),  kinetic  theory 
provides  formulas  for  the  calculation  of  the  transport  properties.  Such 
calculations  have  been  carried  out  by  Amdur  and  Mason  (Ref.  5)  for 
inert  gases  and  mixtures  between  1000  K  and  15,  000  K.  Subsequent 
transport  property  measurements  in  arcs  (Refs.  6,  7)  and  shock  tubes 
(Refs.  8,  9)  have  shown  general  agreement  with  the  calculations  of 
Amdur  and  Mason.  A  second  type  of  practical  information  provided  by 
kinetic  theory  is  the  relationship  between  the  various  measurable  gas 
properties.  A  good  example  is  the  derivation  of  diffusion  coefficients 
of  a  mixture  from  the  viscosity  (Ref.  10)  or  thermal  conductivity 
(Ref.  11)  of  the  mixture  and  the  pure  gases.  Again,  for  the  inert  gases 
the  mixing  rules  from  kinetic  theory  are  in  agreement  with  experiment 
within  a  few  percent. 

The  transport  of  energy  in  polyatomic  gases  must  include  energy 
transported  by  internal  modes  such  as  molecular  rotation.  Recent  ^ 
calculations  of  transport  properties  of  polyatomic  gases  above  2000  K 
proceed  in  two  steps.  First,  using  the  kinetic  theory  appropriate  for 
monatomic  gases  and  cross  sections  (from  molecular  beam  experi¬ 
ments,  for  example)  the  transport  due  to  the  translational  motion  of 
the  molecules  is  calculated.  The  effects  of  internal  modes  are  then 
added  based  on  semiempirical  relations  (for  instance,  by  creating  a 
Eucken  factor  for  internal  energies). 
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Recently  the  formal  kinetic  theory  has  been  extended  to  poly¬ 
atomic  gases  through  the  efforts  of  Wang  Chang,  Uhlenbeck  and  DeBoer 
(Ref.  12).  The  resulting  expressions  are  too  complicated  to  serve  as 
a  recipe  for  the  practical  calculation  of  transport  coefficients.  How¬ 
ever,  further  development  due  to  Mason,  Monchick  and  their  associates 
(Refs.  13-16)  has  provided  a  firm  kinetic  theory  basis  for  separating 
the  translational  and  internal  modes  of  transport.  The  most  important 
feature  of  the  kinetic  theory  of  polyatomic  gases  is  the  dependence  of 
energy  transport  on  diffusion  coefficients  and  relaxation  times.  Al¬ 
though  theory  and  experiment  are  not  in  complete  accord  at  the  present 
time  (Ref.  17),  this  information  is  basic  for  the  calculation  of  the  trans¬ 
port  properties  of  polyatomic  gas  mixtures  such  as  air. 

At  sufficiently  high  temperatures  all  gases  become  mixtures  due 
to  chemical  reactions  such  as  dissociation  and  ionization.  The  presence 
of  chemical  reactions  in  the  gas  leads  to  an  increase  in  thermal  conduc¬ 
tivity  (Ref.  18)  which  is  described  as  a  reaction  conductivity.  The 
reaction  conductivity  results  from  dissociated  or  ionized  species  dif¬ 
fusing  to  low  temperature  regions  and  recombining.  The  heat  of  re¬ 
combination  is  released  in  the  low  temperature  region  resulting  in 
energy  transport.  This  effect  may  be  very  large.  For  instance,  in  , 
air  at  7000  K  the  thermal  conductivity  including  reaction  conductivity 
is  a  factor  of  ~10  greater  than  frozen  thermal  conductivity  (Ref.  19) 
(translational  plus  internal). 

Pressure  affects  the  transport  properties  of  chemically  reacting 
gases  both  by  shifting  the  chemical  equilibrium  and  by  causing  multiple 
collision  effects.  Shifting  the  chemical  reaction  changes  the  transport 
properties  by  changing  the  composition  of  the  mixture.  In  addition,  the 
reaction  conductivity  is  affected  by  pressure.  For  instance,  at  7000  K 
the  thermal  conductivity  in  air  decreases  by  a  factor  of  2  (Ref.  19)  when 
pressure  is  increased  to  30  atm.  This  is  due  to  the  reaction  conduc¬ 
tivity.  The  effect  of  multiple  collisions  on  the  transport  properties  at 
the  same  temperature  and  pressure  on  the  other  hand  is  negligible. 

Ionization  introduces  additional  problems  in  the  determination  of 
transport  properties  both  in  theory  and  experiment.  Electrons  are  very 
efficient  at  exciting  internal  modes.  There  is  also  charge  exchange 
between  ions  and  neutrals  (Ref.  20).  Because  the  coulomb  force  is  long 
range  it  is  necessary  to  consider  multiple  collisions  even  at  low  den¬ 
sities.  A  Chapman- Enskog  procedure  has  been  developed  to  handle 
ionized  gases.  However,  higher  Chapman-Enskog  approximations'  may 
be  necessary  to  calculate  satisfactory  transport  properties  (Refs.  3,21). 
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Two  additional  transport  properties  become  important  in  ionized 
gases,' ,'radiative  flux  and  electrical  conductivity.  For  instance,  the 
caScalde  arc  techniques  of  Burhorn  (Ref.  22)  and  Maecker  and  Wienecke 
(R%f.‘  23)  require  both  of  these  properties  as  inputs.  Other  techniques 
for  measuring  transport  properties  at  temperatures  where  ionization  is 
important  are  affected  at  least  by  radiative  transport.  The  electrical 
cohducti'vity  in  some  cases  is  the  only  experimental  check  possible  on 
tfan'sport  property  calculations. 

Finally,  the  effect  of  pressure  (density)  on  transport  properties 
as  well  as  the  thermodynamics  of  ionized  gases  is  open  to  serious 
qiie’stion.  The  problem  arises  because  it  is  difficult  to  define  what 
is  meant  by  a  free  electron,  since  for  sufficiently  high  densities  the 
electron  always  has  a  potential  energy  equivalent  to  its  thermal  energy. 
In  addition,  radiative  transport  becomes  very  important  at  the  higher 
densities. 

The  temperature-pressure  range  of  available  transport  data  is 
limited  by  the  experimental  methods.  The  classical  methods  such  as 
capillary  flow  or  thin  wire,  require  friction  or  heat  transfer  on  a  solid 
surface.  These  methods  are  limited  by  the  destruction  temperature  of 
solids.  Depending  on  the  pressure,  the  classical  methods  are  limited 
by  the  strength  of  these  solids  at  elevated  temperature.  Viscosity  at 
temperatures  less  than  2000  K  is  the  best  known  transport  property. 
There  are  uncertainties  in  the  order  of  5  to  10%  in  thermal  conduc¬ 
tivity  in  the  same  temperature  range.  The  reasons  for  the  uncertainty 
in  thermal  conductivity  are  due  to  difficulties  in  handling  convection, 
radiation  and  gas  wall  energy  exchange  (accommodation)  (Refs.  24,25). 
The  radiation  problem  is  aggravated  by  high  pressures.  Diffusion  co¬ 
efficients  have  the  largest  uncertainties,  10%  at  room  temperature. 
Diffusion  coefficients  at  temperatures  above  room  temperature  have 
become  available  only  recently  with  the  advent  of  the  point  source 
technique  (Refs.  26,  27).  More  recently,  ultrasonic  techniques  have 
extended  diffusion  coefficients  to  5000  K  (Ref.  28).  Very  few  measure¬ 
ments  of  diffusion  coefficients  at  high  pressures  are  available. 

^  As  mentioned  above,  the  kinetic  theory  of  gases  with  internal 
fno’des  involves  diffusion  coefficients,  relaxation  times  and  bulk  vis¬ 
cosity;  When  calculating  transport  properties  from  molecular  po¬ 
tentials  or  extrapolating  rare  gas  properties  to  high  densities,  it  is 
convenient  to  separate  the  internal  modes  from  the  translational  modes 
of  transport.  This  can  be  done  if  one  knows  the  appropriate  diffusion 
coefficients  and  relaxation  times.  There  are  several  experiments 
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which  are  capable  of  measuring  either  translational  thermal  conduc¬ 
tivity  or  relaxation  times.  Some  of  these  experiments  are  free  jet 
expansion  (Ref.  29),  thermal  transpiration  (Ref.  30)j  and  ultrasonic 
absorption  (Ref.  31).  While  in  the  first  approximation  the  effects  of 
internal  modes  are  not  a  function  of  density,  some  high  density  mea¬ 
surements  of  these  properties  would  be  interesting. 

During  the  past  year  three  experiments  were  chosen  to  answer 
some  immediate  needs  in  the  area  of  transport  properties  for  engi¬ 
neering  purposes.  A  heat  transfer  experiment  was  investigated  for 
measuring  the  effects  of  pressure  on  thermal  conductivity.  A  schlie'ren 
experiment  was  set  up.  This  experiment  will  provide  the  thermal  con¬ 
ductivity  of  gases  at  high  temperature  (2000°K  <  T  <  11,  000°K),  and 
moderate  pressures  (p  ~300  atm).  Finally,  ultrasonic  absorption  mea¬ 
surements  were  made  in  helium  at  temperatures  up  to  3000  K  and 
pressures  of  2  atm. 
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SECTION  II 
DATA  REVIEW 


2.1  DEVIATION  CHARTS 

A  review  of  available  transport  data  for  He,  Ar,  and  Air 

was  undertaken  to  determine  what  additional  measurements  would  best 
contribute  to  present  knowledge  and  to  assess  the  accuracy  of  extrap- 
olati9n  procedures  for  present  needs.  In  addition,  a  limited  comparison 
was,  made  between  theory  and  experiment  for  mixtures  of  these  gases. 

The  viscosity  and  thermal  conductivity  of  each  of  the  gases  under  investi¬ 
gation  are  presented  in  the  form  of  deviation  charts.  The  temperature 
deviation  charts  are  constructed  by  taking  a  reference  viscosity  from 
available  compilations  and  plotting 


(1) 


vs  temperature,  where  rje  measured  viscosity  and  is  the  ref 

erence  viscosity.  The  excess  viscosity,  charts  are  plots  of 


rj 


+ 


(2) 


vs  pressure  where  tJq  is  the  low  pressure  limit  (measured  in  the  same 
experiment  where  possible).  The  t|q' s  are  included  in  the  temperature 
deviation  chart  for  comparison  with  other  low  pressure  measurements. 
The  pressure  was  used  instead  of  the  more  meaningful  variable  density 
(Ref.  32)  because  in  applications  pressure  is  given  more  often  than 
density.  The  notation  for  the  other  transport  coefficients  is  analogous. 

The  deviation  chart  presentation  is  chosen  because  it  allows  rapid 
determination  of  the  accuracy  of  theory  relative  to  experiment.  The 
consistency  of  measurements  from  different  sources  may  be  determined 
by  inspection.  In  general,  the  internal  consistency  of  measurements  is 
much  better  than  the  agreement  between  theory  and  experiment.  Finally, 
transport  properties  for  use  in  engineering  applications  may  be  obtained 
by  adding  the  deviations  to  the  reference  values.  For  instance,  for 
viscosity 
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T)  =  rj  (1  +  A  )  .  (3) 

•*^7]  ■  . , 

The  deviation  charts  may  be  continuously  updated  as  new  measurements 
become  available. 


2.2  HELIUM 

The  viscosity  of  helium  has  been  measured  at  pressures  up  to 
820  atm  (Ref.  33)  and  temperatures  up  to  2340°K  (Ref.  34)  by  the  cap¬ 
illary  flow  technique.  Additional  measurements  have  been  made  using 
ultrasonic  absorption  (Ref.  6)  at  temperatures  up  to  ~4000  K  in  the 
experimental  part  of  the  present  effort.  A  deviation  chart  of  the. nor¬ 
malized  viscosity  of  helium  vs  temperature  is  shown  in  Figure  1  from 
several  recent  investigations.  The  reference  viscosity  was  taken  from 
the  work  of  Lick  and  Emmons  (Ref.  35).  The  low  temperature 
(T  <  3000  K)  viscosities  from  the  various  sources  agree  to  within  2  to 
7%.  The  ultrasonic  measurements  have  a  small  uncertainty,  ~2%,  be¬ 
cause  Prandtl  numbers  based  on  the  work  of  Lick  and  Emmons  had  to 
be  used  to  reduce  the  sound  absorption  to  viscosity. 

The  average  of  the  measured  viscosities  is  about  2%  higher  than 
the  theoretical  calculations  of  Lick  and  Emmons.  At  the  high  tempera¬ 
ture  the  viscosities  of  Amdur  and  Mason  (Ref.  5)  deviate  from  the  values 
of  Lick  and  Emmons  by  about  12%.  Recently  revised  scattering  resialts 
published  by  Amdur  and  Jordan  (Ref.  4)  should  give  viscosities  which 
are  even  closer  to  the  values  of  Lick  and  Emmons.  Some  calculations 
by  Devoto  and  Li  (Ref.  3)  have  been  carried  out  above  4000°K  but  the 
resulting  viscosities  are  systematically  3%  higher  than  the  values  of  , 
Lick  and  Emmons  based  on  the  same  high  energy  potentials  (Ref.  36). 

The  interaction  potential  of  Mason  and  Rice  (Ref.  37)  at  low  tempera¬ 
ture,  T  <  2000  K,  is  probably  the  best  description  of  the  viscosity  of 
helium.  At  temperatures  between  2000  K  and  10,  000  K  the  tabulation 
of  Lick  and  Emmons  should  be  used.  For  ionized  helium  the  calcula¬ 
tions  of  Devoto  and  Li  appear  to  be  the  most  reliable. 

The  thermal  conductivity  of  helium  (at  1  atm)  vs  temperature  de¬ 
viation  chart  is  shown  in  Figure  2.  The  thermal  conductivity  has  been 
measured  at  temperatures  up  to  2400  K  (Ref.  38)  using  the  thin  wire 
technique.  Measurements  up  to  this  temperature  have  also  been  made 
by  Blais  and  Mann  (Ref.  39)  with  an  adaptation  of  a  thermal  diffusion 
column.  These  measurements  are  most  likely  too  high  by  about  20% 
at  the  highest  temperature  (Refs.  37,40)  (~2000  K).  Measurements  at 
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temperatures  up  to  about  600°K  and  pressures  up  to  500  atmospheres 
were  recently  made  (Ref.  41).  The  reference  thermal  conductivities 
are  again  taken  from  Lick  and  Emmons. 

The  experimental  measurements  of  thermal  conductivity  before 
1962  have  been  reviewed  by  Ho  and  Leidenfrost  (Ref.  42).  The  experi¬ 
mental  data  is  consistent  within  5%  at  temperatures  up  to  2400  K.  The 
overall  behavior  of  the  thermal  conductivity  is  similar  to  the  viscosity 
with  the  exception  that  the  deviations  amongst  data  from  various  sources 
are'about  5%  greater. 

'  Some  high  temperature  measurements  of  the  temperature  depend¬ 
ence  of  thermal  conductivity  (Ref.  43)  using  the  heat  transfer  to  the  end 
wall  of  a  shock  tube  are  shown  in  Figure  2.  These  measurements  de¬ 
viate  -from  Amdur  and  Mason*  s  (Ref.  5)  thermal  conductivity  data  by  up 
to  30%.  However,  the  thermal  conductivity  from  the  eiid  wall  experi¬ 
ment  agrees  more  closely  with  thermal  conductivity  calculated  from  the 
recent  theoretical  He-He  cross  sections  of  Phillipson  (Refs.  3,  36)  at 
temperatures  above  4000  K.  This  is  of  questionable  significance  since 
the  shock  tube  data  reported  does  not  contain  significant  information 
about  thermal  conductivity  above  4000°K.  In  fact,  the  significance  of 
the  end  wall  data  is  weighted  toward  the  lower  temperatures.  However, 
the  low  temperature  data  are  in  general  agreement  with  the  potential 
functions  of  Mason  and  Rice.  They  were  able  to  correlate  all  the  exist-  . 
ing  data  (up  to  1954)  for  viscosity,  thermal  conductivity  and  second 
virial  coefficients  for  He  to  within  2%. 

The  viscosity  and  thermal  conductivity  of  helium  show  only  a  slight 
pressure  dependence.  For  this  reason  definitive  experimental  measure- 
mehts  of  this  effect  have  only  recently  become  available. 

The  deviation  chart  for  viscosity  vs  pressure  for  several  investi¬ 
gators  is  shown  in  Figure  3.  The  viscosity  increases  only  5%  with  a 
300  atm  pressure  increase  at  100  C.  The  Chapman- Enskog  (Ref.  2) 
theory  gives  results  which  agree  rather  well  with  available  data  (Ref. 

32).  Values  of  A  ,  B  and  D  which  fit  the  equation 
\  \  \ 

j 

2  2 

rj  =  r]  +  A  p  +  B  p  In  p  +  D  p  +  (4) 

I  O  \  A.  \ 


have  been  investigated  by  Sengers. 
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Similar  results  for  the  thermal  conductivity  of  helium  are  shown 
in  Figure  4.  The  thermal  conductivity  at  high  pressures  has  been  re¬ 
viewed  by  Ho  and  Leidenfrost.  They  found  that  the  expression  of  Snider 
and  Curtiss  (Ref,  44)  agreed  with  these  data  much  better  than  the 
Fnskog  theory. 


2.3  ARGON 

The  viscosity  of  argon  at  atmospheric  pressure  has  been  reviewed 
by  several  authors  (Refs.  45,  46,47).  Argon  measured  viscosities  up  to 
1800  K  are  consistent  to  within  1%.  However,  several  efforts  have 
failed  to  find  a  simple  potential  function  which  reproduces  the  viscosity 
within  fee  experimental  uncertainties  over  the  temperature  range  700°K 
to  1800  K.  Extrapolation  of  experimental  viscosities  of  argon  to  4000°K 
using  a  potential  function  based  on  low  temperature  data  leads  to  errors  ■ 
of  the  order  of  50%.  High  temperature  measurements  of  argon  viscosity 
have  been  made  using  the  ultrasonic  technique  at  temperatures  up  to 
7000  K.  At  the  highest  temperature,  the  measured  thermal  conductivity 
from  the  cascade  arc  experiment  of  Knopp  and  Cambel  (Ref.  7)  was  used 
to  obtain  viscosity.  A  deviation  chart  for  argon  viscosity  is  shown  in 
Figure  5.  The  measured  viscosity  above  1000  K  agrees  with  the  Amdur 
and  Mason  (Ref,  5)  calculations  up  to  temperatures  where  ionization' 
occurs,  ~6000  K. 

The  thermal  conductivity  of  argon  has  been  measured  over  a  much 
larger  temperature  range  than  any  other  gas.  The  data  from  most  in¬ 
vestigations  below  3000  K  is  consistent  to  within  3%  as  may  be  seen  by 
inspecting  the  deviation  chart  in  Figure  6.  The  scatter  in  the  data  from 
each  investigation  for  T  <  2000  K  shows  that  thermal  conductivity  data 
is  not  as  self-consistent  as  the  viscosity  measurements. 

The  thermal  conductivity  of  argon  in  the  temperature  range  from 
3000  K  <  T  <  12,  000°K  is  shown  in  Figure  7.  The  thermal  conductivity 
has  been  measured  by  the  cascade  arc  technique  (Refs.  7,48,49,  50)  and 
an  ultrasonic  technique  (Refs.  6,  8,  9).  The  measurements  are  compared 
to  the  calculations  of  Devoto  (Ref.  21)  and  Yos  (Ref.  19).  The  Maecker 
type  measurements  contain  the  translation  thermal  conductivity  and  the 
reaction  thermal  conductivity  (equilibrium  conductivity).  The  average 
thermal  conductivity  falls  about  10%  above  theory  with  a  possible  un¬ 
certainty  of  about  +  30%  due  to  the  lack  of  accurate  electrical  conduc¬ 
tivity  and  radiative  properties  for  argon.  The  ultrasonic  data  measures 
the  translational  (frozen)  thermal  conductivities  of  argon.  The  reference 
thermal  conductivity  for  the  cascade  arc  measurements  is  the  equilibrium 
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value  and  for  the  ultrasonic  method  is  the  translational  value.  The 
largest  contribution  to  the  sound  absorption  in  argon  between’ 7000  and 
12,  000  K  is  the  thermal  conductivity.  Therefore,  the  viscosity  from 
Devoto' s  theory  was  used  with  the  sound  absorption  to  obtain  the  ther¬ 
mal  conductivity  data  points.  The  ultrasonic  measurements  indicate 
that  Devoto' s  viscosity  as  well  as  thermal  conductivity  are  accurate. 

The  viscosity  of  argon  at  elevated  densities  has  been  measured 
at  pressures  up  to  about  2000  atm.  The  excess  viscosity  chart  for 
argon  vs  pressure  is  shown  in  Figure  8.  The  pressure  effect  in  argon 
is  much  larger  than  helium,  about  a  15%  increase  with  a  pressure  in- 
crease  of  100  atm  at  25  C.  The  thermal  conductivity  of  argon  increases 
somewhat  faster  than  viscosity,  30%  in  100  atm.  An  excess  thermal 
conductivity  plot  is  drawn  up  in  Figure  9.  The  transport  properties  of 
high  pressure  argon  have  been  compared  to  a  modified  Chapman- Enskog 
theory. , (Ref.  33^  and  are  within  12%  of  each  other  over  the  temperature 
range'  from  250  K  to  375°K  and  from  100  to  700  atm. 

2.4  POLYATOMIC  GASES 

■Hie  review  of  the  polyatomic  gas  data  requires  some  new  con¬ 
cepts  largely  to  account  for  the  transport  of  energy  by  internal  modes. 
Similar  considerations  must  be  applied  to  monatomic  gases  at  high 
temperatures  to  describe  energy  transport  due  to  exchange  between 
translation  and  electronic  states.  The  viscosity  depends  on  the  internal 
mode  only  to  the  extent  that  the  cross  sections  depend  on  the  internal 
state  of  the  particle.  Thus,  for  the  calculation  of  viscosity  the  gas  is  a 
mix^re  of  different  particles  each  corresponding  to  a  degree  of  excita¬ 
tion^  ‘The  thermal  conductivity  on  the  other  hand  may  be  expressed  as 
a  sum  of  a  translational  (analogous  to  viscosity)  and  an  internal  thermal 
conductivity.  The  internal  thermal  conductivity  describes  an  excited 
molecule  diffusing  into  a  cooler  region  of  the  gas  and  giving  up  its  in¬ 
ternal  energy  to  the  neighboring  cooler  molecules  at  a  rate  depending 
on  the  relaxation  time  for  the  exchange  process  (Ref.  17). 

The  problem  of  interncil  modes  (Ref.  13)  was  first  approached  by 
considering  deviations  of  the  Eucken  factor  (units  chosen  to  make  f 
dimensionless) 


(5) 
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from  the  theoretical  value  for  monatomic  gases.  The  Eucken  factor  for 
monatomic  gases  as  a  function  of  temperature  (Ref.  26)  is  shown  in 
Figure  10.  An  inverse  power  law  interaction  potential  gives  a  Eucken 
factor  which  varies  from  2.  522  (n  =  5)  to  2.  50  (n  =  oo  for  a  hard  sphere). 
The  addition  of  an  attractive  potential  to  the  hard  sphere  interaction 
causes  the  Eucken  factor  to  vary  with  temperature.  Examination  of 
Figure  10  shows  monatomic  gases  have  Eucken  factors  which  fall  within 
8%  of  the  simplest  theory.  In  addition,  there  are  experimental  indica¬ 
tions  that  the  Eucken  factor  systematically  increases  with  temperature 
and  atomic  number  for  the  monatomic  gases  (Ref.  25).  This  should  be 
expected  if  f  depends  only  on  the  interaction  potential.  The  Eucken 
factor  for  polyatomic  gases  on  the  other  hand  deviates  typically  by  50% 
from  the  classical  value  of  f.  The  f  for  polyatomic  gases  depends  in  a 
very  complicated  way  on  temperature  and  molecular  properties  (Refs. 

13,  51).  In  general,  the  Eucken. factor  for  polyatomic  gases  is  expressed 
as  the  sum  of  an  internal  and  a  translational  part  so  that 


\ 

0 


+  C 


f.  ^ 

V.  .  int 
int 


(6) 


The  value  of  the  translational  Eucken  factor,  f|.y,  may  be  obtained  from 
viscosity  or  diffusion  data.  The  value  of  the  internal  Eucken  factor  on 
the  other  hand  must  be  evaluated  from  theory  unless  both  \  and  rj  are 
known. 

As  mentioned  in  the  Introduction,  the  exact  kinetic  theory  for 
polyatomic  gases  has  been  developed  in  recent  years.  The  resulting 
relations  will  be  used  below  to  analyze  experimental  transport  data  in 
hydrogen  and  nitrogen.  The  expressions  obtained  directly  from  the 
Boltzmann  equation  for  polyatomic  gases  (Refs.  13,  15)  are  used  below. 
For  a  gas  with  a  single  internal  degree  of  freedom  the  internal  thermal 
conductivity  is 
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and  is  the  diffusion  coefficient  for  the  internal  mode,  is  the 

specific  heat  for  that  mode,  and  Z  is  the  average  number  of  collisions 
necessary  to  exchange  internal  and  translational  energy.  In  most  cases 
the  denominator  of  F  may  be  neglected.  The  translational  thermal  con¬ 
ductivity  is 


The  Eucken  factor  may  be  found  from  (7)  and  (9)  if  the  collision  number, 
viscosity  and  internal  diffusion  coefficient  are  known. 

*  -  I 

■^In  practice  rotation  and  electronic  excitation  are  the  only  impor¬ 
tant  degrees  of  freedom  for  the  calculation  of  thermal  conductivity  from 
(7)  and  (9).  Since  accurate  experimental  data  is  available  only  up  to 
about  2000  K  only  rotation  need  be  considered.  The  internal  diffusion 
coefficient  for  rotation  may  be  taken  to  be  the  self -diffusion  coefficient. 


For  most  gases  the  quantity 


pD 


11  . 


is  ~1.  13. 


The  rotational  col¬ 


lision  numbers  may  be  calculated  using  an  expression  due  to  Parker 
(Ref,  52)  with  constants  adjusted  to  reproduce  experimental  at 

room  temperature.  Theoretical  Eucken  factors  based  on  the  above 
assumptions  agree  with  the  experiment  for  most  gases  within  the  ex¬ 
perimental  error  (Refs.  13,  51), 


■  These  concepts  may  now  be  applied  in  the  analysis  of  the  experi¬ 
mental  data  for  H2  and  N2  . 


2.5.  HYDROGEN 

The  transport  properties  of  hydrogen  are  not  as  well  known  as  the 
other  common  gases.  '  The  deviation  chart  for  the  viscosity  of  H2  vs 
temperature  at  one  atmosphere  is  shown  in  Figure  11.  During  the  past 
few  years  the  viscosity  of  hydrogen  has  been  measured  up  to  2400  K. 
The  early  results  are  consistent  within  2%.  The  most  recent  measure¬ 
ments  (Ref.  34)  are  systematically  6%  higher  than  earlier  work  of 
Trautz  et  al.  (Refs.  53,  54)  and  the  tabulation  of  Hilsenrath  (Ref.  55). 

At  temperatures  above  1300  K  the  viscosities  of  Guevara  and  Wageman 
are  systematically  higher  than  the  viscosities  based  on  the  potential 
function  of  Mason  and  Rice  by  about  10%.  Heat  transfer  measurements 
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on  the  end  wall  of  the  shock  tube  are  suited  to  resolve  the  above  men¬ 
tioned  discrepancy. 

A  very  limited  number  of  viscosity  measurements  for  hydrogen 
are  available  at  high  pressures,  the  most  recent  being  the  Russian' 

(Ref.  56)  measurements  relative  to  nitrogen  at  temperatures  up  to 
1000  K  and  pressures  up  to  about  500  atm.  Sengers  (Ref.  32)  has  com¬ 
pared  the  viscosity  of  hydrogen  up  to  800  amagats  to  the  Enskog  theory. 

The  thermal  conductivity  of  hydrogen  has  been  studied  extensively 
at  temperatures  up  to  2400  K  at  one  atmosphere.  Thermal  conductivity 
measurements  have  been  attempted  up  to  3000  K  by  a  new  technique 
(Ref.  57).  However,  further  experiments  have  to  be  done  to  assess  the 
errors  in  this  system. 

The  deviation  chart  for  thermal  conductivity  of  hydrogen  vs  tem¬ 
perature  (Figure  12)  shows  a  large  amount  of  scatter,  6%,  in' the  mea¬ 
surements  from  different  investigators  at  temperatures  up  to  1000  K. 
Higher  temperature  data  shows  even  less  consistency.  The  data  of 
Blais  and  Mann  are  too  high  by  about  20%  at  the  highest  temperatures. 
This  same  discrepancy  was  pointed  out  in  the  helium  measurements  of 
Blais  and  Mann  above. 

t 

The  translation  thermal  conducti'vity  is  somewhat  difficult  to.  cal¬ 
culate  for  hydrogen.  The  Eucken  factor  for  hydrogen  (Figure  13)  does 
not  behave  as  expected  at  high  temperatures,  Saxena  (Ref,  40)  and 
Mason  (Ref.  37)  conjectured  that  the  difficulties  were  due  to  errors  in 
the  thermal  conductivity  measurements.  The  most  recent  measure¬ 
ments  seem  to  bear  this  out.  Fortunately,  the  internal  thermal  con¬ 
ductivity  of  hydrogen  is  a  much  smaller  effect  than  it  is  in  other  diatomic 
gases. 

2.6  NITROGEN 

The  transport  properties  of  N2  are  known  better  than  any  other 
gas  with  the  possible  exception  of  Ar.  The  viscosity  has  been  measured 
(Ref.  58)  up  to  1400  K  with  a  majority  of  the  points  within  1%.  Near 
room  temperature  several  investigators  have  obtained  the  same  results 
to  within  a  tenth  of  one  percent.  The  deviation  chart  for  these  measure¬ 
ments  is  shown  in  Figure  14.  Some  measurements  of  viscosity  by 
Wienecke  in  an  electric  arc  are  available.  However,  these  values  are 
very  different  from  theoretical  results.  Olher  measurements  must  be 
done  to  confirm  the  work  of  Maecker  and  Wienecke  (Ref.  23), 
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The  viscosity  of  nitrogen  at  high  pressure  has  been  investigated 
up  to  820  atmospheres  at  temperatures  up  to  500  K  (Ref.  33).  The 
pressure  effect  is  quite  large,'  the  viscosity  almost  doubling  in  10^  at¬ 
mospheres.  The  deviation  chart  for  the  available  high  density  nitrogen 
viscosity  was  compared  to  the  modified  Chapman- E ns kog  (Ref.  32) 
relation  and  found  to  be  within  4%.  The  available  excess  viscosity  data 
are.  shown  in  Figure  15. 

The  thermal  conductivity  of  nitrogen  has  been  measured  up  to 
1400  K  at  atmospheric  pressure.  The  earliest  work  tabulated  by 
Hilsenrath  et  al.  showed  a  self-consistency  within  1%  up  to  600  K.  At 
higher  temperatures  more  recent  data  showed  scatter  of  10%.  Bargaftik 
and  Zimina  (Ref.  58)  reviewed  the  available  thermal  conductivity  data  at 
high  temperatures,  applying  corrections  for  temperature  jump  (Ref.  25) 
to  the  measurements  of  several  authors,  and  they  found  the  available 
data  consistent  to  within  2%.  These  corrections  reduced  the  scatter  in 
the  Eucken  factor  vs  temperature  curves  for  N2  and  brought  the  theory 
of  .Mason  et  al.  into  much  closer  agreement  with  experiment.  A  devia¬ 
tion  chart’ for  thermal  conductivity  is  shown  in  Figure  I6. 

..  The  Eucken  factor  vs  temperature  plot  is  given  in  Figure  17.  The 
rotation  numbers  up  to  1500  K  for  this  graph  were  taken  from  jet  expan¬ 
sion  experiments  of  Miller  and  Andres  (Ref.  .29)  and  ultrasonic  measure¬ 
ments  (Ref.  31).  The  temperature  variation  of  the  Eucken  factor  of 
nitrogen  is  an  excellent  test  of  the  polyatomic  gas  theory  because  the 
necessary  numbers  are  available. 

.  The  thermal  conductivity  of  nitrogen  has  been  measured  in  the 
constricted  arc  at  one  atmosphere  and  temperatures  above  5000  K  by 
several  investigators  (Refs.  23,  50,  59i  60).  A  deviation  chart  of  ex¬ 
perimental  high  temperature  thermal  conductivity  values  vs  temperature 
is  shown  in  Figure  18.  The  recent  calculations  of  Yos  (Ref.  19)  are  used 
as  the  reference  values.  The  original  data  of  Maecker  (Ref.  24)  are  in 
good  agreement  with  theory  between  5000  and  6000  K.  However,  .from 
7000  to  15,  000  K  the  measured  values  are  high  and  depend  on  arc  current. 
The  recent  work  at  Avco  (Ref.  6I)  showed  that  a  reasonable  radiative  cor¬ 
rection  broughtthe  Maecker  data  in  reasonable  agreement  with  theory,  and 
eliminated  the  current  dependence.  The  same  correction  applied  to  the 
thermal  conductivity  measurements  at  Avco  also  eliminates  the  current 
dependence.  However,  the  Avco  data  are  about  50%  above  the  Maecker 
data  and  the  calculations  of  Yos. 

The  theoretical  uncertainties  in  these  measurements  are  largely 
due  to  uncertainties  in  the  calculated  radiation  and  electrical  properties 
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of  gases  which  are  necessary  to  reduce  the  data.  Again  it  should  be 
noted  that  near  7000  K  and  again  near  15,  000°K  the  reaction  thermal 
conductivity  dominates  over  the  translational  thermal  conductivity. 

Only  at  about  10,  000  K  do  these  experiments  have  any  information  on 
the  translational  thermal  conductivity.  The  errors  associated  with  the 
cascade  arc  technique  are  further  discussed  below  in  the  methods  section. 


2.7  AIR 


The  viscosity  of  air  has  been  the  subject  of  a  recent  review  by 
Hing  Yi  Lo  et  al,  (Ref.  62).  The  experimental  results  for  viscosity 
vs  temperature  are  very  consistent,  2%  up  to  temperatures  of  1000  K 
and  1%  up  to  2000  K.  This  is  shown  in  the  deviation  chart  in  Figure  19. 
The  reference  viscosity  for  the  chart  is  taken  from  the  correlation  for¬ 
mula  given  in  Ref.  62.  The  pressure  data  correlates  especially  well  at 
high  pressures,  while  the  scatter  at  lower  pressures  is  mainly  due  to 
lack  of  significance  in  the  excess  viscosity.  An  excess  viscosity  chart 
is  drawn  up  for  viscosity  vs  pressure  in  Figure  20.  The  correlation 
formula  from  Ref.  62  was  used  for  the  low  density  viscosities. 

The  mixing  rules  for  the  thermal  conductivity  of  a  polyatomic  gas 
mixture  such  as  air  are  not  fully  understood.  Part  of  the  difficulty  is 
that  insufficient  information  on  collision  numbers  and  diffusion  coeffi¬ 
cients  (especially  the  latter)  exists  to  test  the  relative  merits  of  various 
theories.  Another  problem  is  the  complexity  of  the  exact  theory.  The 
available  experimental  data  for  mixtures  involving  polyatomic  gases 
have  been  compared  to  theory  by  several  authors  (Refs.  16,  63).  The 
expression  due  to  Hirschfelder  (Ref.  64),  seems  to  agree  with  experi¬ 
ment  better  than  the  recent  attempts  to  treat  polyatomic  mixtures  more 
exactly  (within  10%  in  the  worst  case).  Hirschfelder' s  expression  is 


\  . 
mix 


1  + 


JL  D..  X. 

/,  D..X. 

^  ij  1 

j  =  1 


(10) 


where  \  is  the  thermal  conductivity  of  the  mixture  with  frozen 

mix  o 

internal  modes  (translational  thermal  conductivity).  and  \  ,  are 

the  total  and  translational  thermal  conductivity,  respectively,  of  the 
pure  gases.  Djj  and  are  the  self-diffusion  and  binary  diffusion  co¬ 
efficients,  respectively.  The  X^  are  the  mole  fractions.  The  collision 
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numbers  enter  into  \  mainly  through  \  . 

'  mix  ■  1 

The  most  appropriate  way  to  calculate  the  thermal  conductivity  of 
a  mixture  is  to  obtain  the  thermal  conductivity  of  the  pure  gases  exactly 
from  Eq.  (9),  and  to  use  these  values  in  Eq.  (10).  The  inverse  problem 
arises  when  the  total  thermal  conductivity  is  given  and  the  translational 
thermal  conductivity  must  be  found.  For  instance,  if  we  wish  to  apply 
a  density  correction,  the  internal  thermal  conductivities  of  the  pure 
components  \  \  j  may  be  calculated  at  one  atmosphere,  and  the  in¬ 

ternal  thermal  conductivity  of  the  mixture  calculated  from 


It  could  be  noted  that  because  of  the  very  unsatisfactory  state  of  the 
mixing'  rules  at  high  pressures,  internal  thermal  conductivities  should  > 
be  calculated  from  low  density  data. 

1  ■ 

The  thermal  conductivity  of  air  is  not  as  well  known  as  the  vis¬ 
cosity.  The  early  measurements  tabulated  by  Hilsenrath  up  to  about 
600  K  show  about  5%  scatter.  Tsederberg  (Ref.  25)  reports  some  data 
which  he  tabulates  eifter  some  modification  up  to  1300  K.  In  the  same 
reference  thermal  conductivity  measurements  by  other  Russian  workers 
are  quoted  up  to  400  atm.  However,  Tsederberg  noted  that  their  low 
pressure  values  for  nitrogen  differed  by  6-9%  from  other  reliable  values. 
A  deviation  chart  for  the  thermal  conductivity  of  air  is  shown  in  Figure 

21: 

Eucken  factors  have  been  calculated  from  the  data  in  Figures  19 
and  21  and  are  somewhat  lower  than  direct  measurements  of  Eckert 
et  al.  (Ref.  65).  A  deviation  chart  for  the  pressure  effect  for  the  ther¬ 
mal  conductivity  of  air  is  shown  in  Figure  22.  The  experimental  data 
are  compared  to  a  calculated  thermal  conductivity.  The  internal  ther¬ 
mal  conductivity  was  calculated  from  Eq.  (11)  and  the  data  for  pure 
oxygen  and  nitrogen.  The  translational  thermal  conductivity  was  then 
found  by  subtracting  the  interned  contribution  calculated  above  from  the 
measured  thermal  conductivity  at  one  atmosphere.  The  modified 
Chapman- Enskog  theory  was  then  used  to  extrapolate  the  translational  ' 
thermal  conductivity  to  high  pressures.  The  internal  component  was 
added  again  to  the  extrapolated  translational  thermal  conductivity  to 
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obtain  the  total  thermal  conductivity  at  high  pressures.  The  calculated 
and  experimental  values  are  within  16%. 

An  experimental  determination  (Ref.  66)  of  thermal  conductivity 
of  air  between  1000  K  and  6500  K  has  recently  been  made  using  a 
schlieren  interferometer  technique  on  the  end  wall  boundary  layer 
of  the  shock  tube.  The  measurements  were  made  at  a  constant  density 
of  one-tenth  room  temperature  atmospheric  density.  The  pressure 
range  of  the  data  varies  from  one-third  atmosphere  at  1000  K  to  about 
3  atmospheres  at  6000  K,  A  deviation  chart  for  these  measurements 
is  shown  in  Figure  23.  The  reference  thermal  conductivities  were  taken 
from  Hansen  (Ref.  67).  Also  shown  are  the  recent  calculations  of  Yos 
(Ref.  19).  The  quoted  uncertainties  in  the  measurement  are  50%.  How¬ 
ever,  this  is  greater  than  the  maximum  error,  the  average  of  which  is 
probably  a  good  representation  of  the  thermal  conductivity  of  air.  At 
temperatures  below  2000  K  the  measured  values  agree  with  the  recent 
calculations  of  Yos.  At  higher  temperatures  the  experimental  data  varies 
by  up  to  a  factor  of  two  from  the  two  theoretical  values,  depending  on  the 
temperature.  The  uncertainties  at  this  time  in  the  thermal  conductivity 
of  air  above  2000  K  are  as  large  as  +  75%,  especially  near  temperatures 
where  dissociation  and  ionization  are  about  50%  complete.  Since  dis¬ 
sociation  and  ionization  depend  markedly  on  pressure  between  1  and  10^ 
atm,  this  uncertainty  exists  over  the  entire  temperature  range  of  2000  K- 
20,  000  K.  Until  the  pressure  effects  of  chemical  reactions  on  the  ther¬ 
mal  conductivity  of  air  are  better  known,  consideration  of  higher  order 
collisions  will  not  improve  the  situation.  Fortunately,  at  the  higher  tem¬ 
peratures,  T  >  3000  K,  and  at  a  fixed  pressure,  the  density  is  reduced 
by  factors  of  10  to  100  over  the  density  at  room  temperature.  At  the 
higher  temperatures,  pressures  of  10^  to  lO'^  atm  are  needed  to  cause 
up  to  30%  increase  in  thermal  conductivity  due  to  triple  collisions. 
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SECTION  III 
METHODS 


3.1  '  CLASSICAL  TECHNIQUES 

'  The  problem  of  measuring  transport  properties  reduces  to  creat¬ 
ing  'flows  which  are  amendable  to  analysis.  This  means  flows  with  high 
degrees  of  symmetry  and  small  variations  in  velocity,  temperature, 
pressure  or  species  concentration.  At  low  temperatures,  T  <  3500°K, 
it' is  feasible  to  confine  the  gas  to  very  small  regions  by  solid  bounda- 
riesl  At  high  temperatures,  T  >  3500  K,  it  is  no  longer  possible  to 
contain  the  gases  with  any  known  material.  The  hot  gases  must  be 
sensed  through  boundary  layers  using  light  and  sound  waves.  Indirect 
information  on  transport  properties  of  a  high  temperature  gas  may  be 
o'btaine'd  by  measuring  the  heat  flux  through  the  boundary  layer.  How¬ 
ever,  even  this  technique  requires  accurate  knowledge  of  the  tenipera- 
ture  and  pressure  of  the  high  temperature  gas.  These  thermodynamic 
quantities  must  then  be  sensed  through  a  thermal  boundary  layer.  The 
bigges't  problem  is  to  obtain  a  uniform,  stable  high  teniperature  gas 
soiirc'e'.  This  problem  becomes  more  severe  as  the  pressure  is  in¬ 
creased.  Development  of  high  temperature  gas  sources  has  proved  to 
be  ^e  most  expensive  and  time  cons\uning  task  in  obtaining  transport 
properties  above  3500  K. 

■  The  classical  techniques  for  measuring  viscosity  and  thermal  con¬ 
ductivity  are  discussed  and  appropriate  references  given  by  Sengers 
(Ref.  68).  The  present  discussion  will  be  limited  to  methods  which  are 
most  suitable  for  extension  to  the  extremes  of  temperature  and  pressure 
(p  <’  10^  atm,  T  ^  3500°K).  The  main  reason  for  further,  work  below 
3500  K  is  to  support  theoretical  developments  in  dense  gas  transport 
properties.  A  secondary  purpose  is  to  obtain  adequate  potential  functions 
for  extrapolation  of  gas  properties  to  higher  temperatures.  Past  experi¬ 
ence  has  shown  that  both  types  of  activities  require  extremely  accurate 
(0.  1%)  measurements  to  be  of  any  use. 

The  capillary  flow  method  is  the  simplest  technique  for  the  mea¬ 
surement  of  viscosity.  The  basic  technique  consists  of  forcing  a  steady 
state  flow  rate,  Q,  through  a  capillary  and  measuring  the  resulting  pres¬ 
sure  drop,  p.  The  essentials  of  a  high  pressure-high  temperature  capil¬ 
lary  viscometer  (in  cutaway  view)  are  shown  in  Figure  24.  The  capillary 
is  a  hole  drilled  in  the  wall  of  the  interior  cylindrical  heater  (Guevara 
and  Wageman  Ref.  34).  A  second  guard  heater  with  a  much  thinner  wall 
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is  used  to  maintain  the  uniform  temperature  of  the  inner  heater.  The 
ends  are  attached  to  water-cooled  copper  flanges  which  are  maintained 
at  room  temperatures.  The  heaters  must  be  sufficiently  long  to  bring 
the  test  gas  up  to  the  wall  temperature  before  it  enters  the  capillary. 

For  high  pressures  the  cell  must  be  surrounded  by  a  cooled  high  pres¬ 
sure  support. 

The  advantage  of  drilling  the  capillary  in  the  wall  of  the  inner 
heater  is  that  uniform  temperatures  are  easily  obtained.  In  addition, 
for  high  pressure  applications  the  geometry  in  Figure  24  has  high 
strength.  The  disadvantage  is  that  the  short  length  of  the  capillary 
requires  large  corrections  to  the  simple  Poiseuille  law.  The  wall 
thickness  of  the  inner  heater  is  limited  by  the  power  necessary  to 
maintain  it  at  a  high  temperature.  A  quarter  inch  wall  thickness  on  a 
6"  long  cylinder  requires  in  the  order  of  tens  of  kilowatts.  Some  new 
techniques  have -been  developed  for  manufacturing  small  bore  refrac¬ 
tory  metals  (Mo,  W)  tubes  by  spraying  the  refractory  metal  on  a  form 
and  etching  the  form  away.  However,  the  short  capillary  with  appro¬ 
priate  corrections  is  the  most  promising  technique  for  obtaining  accurate 
viscosity  measurements  at  temperatures  from  1500  C  to  3000  C  and 
pressures  up  to  10^  atm. 

The  capillary  flow  technique  requires  a  method  for  measuring  the 
pressure  drop  across  the  capillary  and  a  pvimp  to  maintain  a  constant 
flow  rate  and  average  pressure.  The  pressure  drop  of  about  7  cm  Hg 
must  be  measured  by  a  device  which  is  capable  of  withstanding  the  high 
mean  pressures.  The  pvimp  usually  has  a  working  fluid  other  than  the 
gas  which  has  a  low  compressibility  thermal  expansion  and  vapor  pres¬ 
sure.  Mercury  has  become  the  standard  choice.  Some  investigators 
have  used  bellows  to  act  directly  on  the  test  gas  as  a  pump. 


High  pressure  capillary  viscometers  have  recently  been  fully  dis¬ 
cussed  by  Giddings,  Kao  and  Kobayashi  (Ref.  69).  The  generalized 
relation  between  viscosity  and  the  experimental  parameters  is 
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(12) 


The  first  terni  in  brackets  is  essentially  the  simple  Poiseuille  law  with 
some  modifications.  The  radius  of  the  capillary  is  replaced  with  an 
average  radius  over  the  length  of  the  tube 
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(13) 


Since  r]  depends  on  a  most  investigators  go  to  great  pains  to  determine 
the  radius  accurately.  The  length  of  the  capillary  is  increased  by  a  con¬ 
stant  multiple  of  the  average  radius  (  n  a)  to  account  for  the  forces  nec¬ 
essary  to  change  the  velocity  distribution  of  the  incoming  fluid  to  the 
parabolic  distribution  characteristic  of  steady  state  Poiseuille  flow.  The 
correction  for  the  compressibility  of  the  gas  is  expressed  as  the  ratio  of 
the  average  density  P  defined  to  the  density  at  the  entrance  of  the  cap¬ 
illary  p  j.  The  average  density  is  given  by 


(14) 


(15) 


This  correction  is  important  near  phase  changes.  A  correction  factor 
for  the  uniformity  of  the  capillary  |  is  included  which  must  be  deter¬ 
mined  by  calibration  with  a  known  gas. 


The  second  term  in  Eq.  (12)  (the  Hagenbach  factor)  corrects  for 
the  acceleration  the  fluid  experiences  going  from  reservoir  flow  into  the 
capillary  flow.  Further  corrections  for  thermal  expansion  and  compres¬ 
sibility  of  the  capillary  walls  are  not  shown.  Also  an  additive  correction 
for  slip  roughly  equal  to  the  ratio  of  mean  free  path  to  capillary  diameter 
is  important  at  pressures  of  one  atmosphere  and  less. 


Thermal  conductivity  has  proved  to  be  a  much  more  difficult  mea¬ 
surement  to  make  atithe  extremes  of  temperature  and  pressure  than 
viscosity.  Even  aftejr  the  very  difficult  problems  of  obtaining  a  uniform 
region  of  test  gas  anji  correcting  for  end  effects  and  edge  effects  have 
been  accomplished,  Iserious  problems  remain  (Ref.  25).  The  problems 
are  convection,  temperature  jump,  radiation  and  temperature  measure¬ 
ment. 
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Convection  occurs  when  dense  fluids  are  generated  above  lighter 
fluids.  Gravitational  forces  cause  the  dense  fluid  to  fall  and  the  lighter 
gas  to  rise.  The  lighter  fluid  carries  its  heat  with  it  adding  to  the  total 
heat  fl\ix.  The  magnitude  of  the  convective  effects  depends  on  a  dimen¬ 
sionless-quantity,  the  Rayleigh  number  (Refs.  69,70), 

2  3  CJ 

Ra=  Gr*Pr  =  P- 

n 


where  g  is  the  acceleration  of  gravity,  Gr  the  Grassoff  number,  and  Pr 
the  Prandtl  nvimber  (these  are  indicated  by  the  product).  The  new  quan¬ 
tities  in  Eq.  (16)  are  the  acceleration  of  gravity  g,  the  thermal  expansion 
coefficient  p,  the  temperature  difference  AT  and  the  thickness  of  the  gas 
layer  f  .  The  heat  transfer  due  to  convection  is  negligible  for  Rayleigh 
numbers  less  than  500.  For  ordinary  gases,  at  atmospheric  densities 
and  with  typical  temperature  differences,  this  criterion  demands  gas 
paths  of  1-0.  1  mm  in  a  thermal  conductivity  apparatus. 

The  technique  which  is  least  sensitive  to  convection  is  a  horizontal 
parallel  plate  apparatus  with  the  high  temperature  plate  on  top.  This 
technique  has  been  used  for  measurements  near  the  critical  density  where 
convection  is  a  very  large  effect  (Ref.  71).  However,  it  is  difficult  to 
achieve  a  simple  temperature  field  between  parallel  plates.  The  other 
extreme  is  the  technique  of  Blais  and  Mann  illustrated  in  Figure  25.  The 
convection  is  driven  by  maintaining  room  temperature  on  the  outer  cyl¬ 
inder  (a  wire  held  at  ~10^  °K).  The  only  regions  with  significant  tem¬ 
perature  drops  in  the  z  direction  are  near  the  ends  of  the  tube.  Therefore, 
near  the  middle  of  the  diffusion  column  (A  B  in  Figure  25)  the  heat  flux  due 
to  convection. 


V  p  C 

V 


a  T 
a  z 


(17) 


is  zero  and  the  heat  conduction  equation  reduces  to  the  radial  term  only. 
The  well-known  solution  to  the  radial  heat  equation  with  large  temperature 
ranges  (\  if  constant)  is 


Q 


\  (T)  dT, 


(18) 
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where  the  symbols  are  provided  in  Figure  25.  The  power  dissipated,  Q, 
in  the  length  of  wire,  L,  is  measured  over  an  entire  range  of  central 
wire  temperatures,  Tj^.  The  resulting  curve  of  Q  vs  Tj^  may  be  differ¬ 
entiated  to  give,  according  to  Eq.  (18),  the  thermal  conductivity. 


(Tj) 


r 


Zn  L 


(19) 


The  most  interesting  aspect  of  this  experiment  is  that  while  Q  is  sig- 

■  '  dQ 

nificantly  sensitive  to  convection,  is  not.  This  technique  has  been 

suggested  as  a  method  to  extend  thermal  conductivity  measurements  to 
temperatures  from  2000  C  to  the  destruction  temperature  of  tungsten. 

The  earliest  measurements  in  helium  (Ref.  39)  shows  serious  discrep¬ 
ancies  with  other  measurements  and  theory.  However,  Saxena  seems 
to  have  successfully  applied  this  technique  (Refs.  72,  73,  74),  and  is 
currently  extending  it. 

A  more  conventional  technique  for  eliminating  convection  is  to 
minimize  spacing  between  the  cold  and  the  hot  elements,  the  tempera¬ 
ture  difference,  AT,  and  operate  at  low  densities.  Even  the  technique 
of  Blais  and  Mann  is  restricted  to  low  densities  because  at  high  Rayleigh 
numbers  convection  is  unstable  and  turbulence  develops.  However,  low 
pressures  lead  to  a  phenomenon  called  temperature  jump.  Energy  trans¬ 
fer  upon  collision  with  a  wall  is  different  from  the  energy  transfer  of  a 
similar  collision  with  a  gas  molecule.  Therefore,  there  is  a  region  up 
to  30  mean  free  paths  long  where  the  temperature  gradient  is  affected  by 
the  efficiency  of  energy  transfer  with  the  wall.  As  a  result  the  hot  wall 
temperature  is  hotter  than  expected  by  extrapolating  the  gradient  in  the 
gas  to  the  wall.  Similarly,  the  cold  wall  is  colder  than  one  would  expect. 

■  f 

The  measured  temperature  difference  between  the  walls  of  the 
thermal  conductivity  cell  must  be  corrected  to  obtain  the  temperature 
drop,  in  the  gas.  The  corrections  for  the  temperature  jump  in  differ¬ 
ent-types  of  thermal  conductivity  cells  have  been  discussed  by  several 
authors  (Refs,  58,  39).  The  basic  technique  is  to  measure  the  thermal 
conductivity  (or  its  equivalent)  at  different  pressures  and  to  extrapolate 
to  infinite  pressure.  One  has^rto  be  careful  to  avoid  the  curvature  which 
tends  to  develop  at  the  higher  pressures  in  polyatomic  gases.  These 
corrections  are  most  important  in  the  thin  wire  variation  of  the  con¬ 
centric  cylinder  method.  Some  authors  (Refs.  58,  75)  have  shown  the 
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temperature  jvimp  to  account  for  2-7%  errors  in  reported  thermal  con¬ 
ductivity  measurements. 


The  major  problem  at  higher  temperatures  (1000°K  <  T  <  3000°K) 
in  conventional  thermal  conductivity  techniques  is  radiation.  The  radia¬ 
tive  transfer  may  be  noinimized  by  making  the  surface  area  of  the  ther¬ 
mal  conductivity  cell  and  the  temperature  difference,  AT,  as  small  as 
possible.  Mainly  for  this  reason  the  accurate  measurements  of  thermal 
conductivity  at  higher  temperatures  have  been  made  with  the  thin  wire 
technique.  The  diagram  in  Figure  25  shows  the  essential  features  of 
the  thin  wire  apparatus.  The  temperature  difference  between  the  central 
wire  and  the  outer  wall  is  of  the  order  of  10  K  in  typical  systems.  Equa¬ 
tion  (18)  relates  the  thermal  conductivity  to  the  measured  heat  flux 
(electrical  power  dissipated  in  the  length  L)  and  temperature  difference. 
Because  the  temperature  difference  is  small,  \  is  essentially  constant. 
The  integral  in  (18)  reduces  to  the  product  of  the  average  thermal  con¬ 
ductivity  and  the  temperature  difference. 

More  recently  a  variable  path  parallel  plate  thermal  conductivity 
cell  has  been  developed  (Ref.  76)  and  used  to  determine  the  thermal 
conductivity  of  helium.  This  type  of.  apparatus  has  the  advantage  in 
that  extrapolation  of  the  effective  thermal  conductivity  to  zero  plate 
separation  may  be  used  to  eliminate  the  effects  of  radiation.  In  addition, 
as  mentioned  above,  this  configuration  minimizes  convection.  Therefore, 
the  variable  path  parallel  plate  system  holds  promise  for  measurements 
at  higher  temperatures  and  densities. 

Finally,  a  problem  which  arises  in  any  thermal  conductivity  cell 
is  the  accurate  measurement  of  temperature.  The  difference  between 
two  measured  temperatures  is  needed  to  determine  the  thermal  conduc¬ 
tivity.  When  temperature  differences  of  the  order  of  10  K  are  involved 
the  temperatures  must  be  accurate  to  0.  1  K  to  provide  1%  accuracy. 
Thermocotgiles  have  been  used  but  these  are  limited  to  temperatures 
below  2000  C.  The  resistance  of  the  central  wire  provided  accurate 
measurement  of  temperatures  for  most  of  the  measurements  with  the 
thin  wire  technique.  The  transit  time  of  an  ultrasonic  wave  has  been 
used  to  measure  temperatures  up  to  2800  K  in  a  rhenivim  wire  (Ref.  77). 
The  ultrasonic  technique  in  solids  along  with  an  ultrasonic  helium  gas 
thermometer  for  calibration  (Ref.  78)  will  do  much  to  improve  the  con¬ 
ventional  techniques  of  thermal  conductivity  determination. 
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3.2  CASCADE  (MAECKER)  ARC 

The  flow  in  a  low  velocity  arc  discharge  in  a  circular  channel 
(constricted  arc,  Figure  26)  which  is  long  in  comparison  to  its  radius, 
is  strongly  dependent  on  the  thermal  conductivity  of  the  gas.  The  mea¬ 
surement  of  the  flow  in  constricted  arcs  has  been  used  to  measure  the 
thermal  conductivity  of  nitrogen  (Refs.  23,  50,  59,  60,  6l)  and  argon 
(Refs,.  7,48,49).  Some  investigators  have  used  previously  reported 
measurements  of  arc  characteristics  to  deduce  the  thermal  conductivity 
of  nitrogen  (Ref.  79)  and  hydrogen  (Refs.  50,  6l).  Measurements  to  date 
have  been  restricted  to  pressures  of  one  atmosphere  and  temperatures 
between  5000°K  to  about  15,  000°K  in  the  arc. 

In  the  original  technique  of  Burhorn  (Ref.  22)  and  Maecker 
(Ref.  23),  the  gas  temperature  was  measured  spectroscopically  and 
estimates  of  radiation  losses  and  electrical  conductivity  were  made. 

The  estimates  were  checked  by  calculating  the  integrated  radiation  loss 
arid  electrical  conductivity  and  comparing  them  with  measured  values. 
Recent  work  has  been  concerned  with  a  better  determination  of  elec¬ 
trical  conductivity  and  radiative  transport.  The  major  problem  involved 
in  extending  the  arc  experiment  to  higher  pressures  rests  in  determining 
the  contribution  of  radiation  to  the  energy  balance. 

Since  the  constricted  arc  is  usually  operated  with  very  low  flow 
velocities,  the  only  important  transport  process  in  this  case  is  the 
transfer  of  heat  between  the  electrically  heated  gas  in  the  center  of  the 
arc’ and  the  water  cooled  walls  by  thermal  conduction.  The  arc  is  usually 
cylindrically  symmetrical  for  large  enough  length  over  diameter  ratio 
(Li/D  ~15).  When  the  constricted  arc  is  made  up  of  alternate  copper  and 
insulator  discs  as  shown  in  Figure  26,  the  thickness  of  the  insulating 
material  must  be  about  l/20  of  the  thickness  of  the  copper  discs  to  pre¬ 
serve  cylindrical  symmetry.  The  energy  balance  for  the  constricted- 
arc'  is  then  (Ref.  50) 


1  _d_ 
r  dr 


+  <te 


Q 
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(20) 


where  r  is  the  radial  distance  from  the  center  of  the  arc,  T  is  the  tem¬ 
perature,  O  the  electrical  conductivity,  £  the  field  strength  and  the 
radiative  energy  flxuc.  Assxuning  that  O’  and  Q^.  are  known,  and  that  the 
electric  field  intensity  and  temperature  as  a  function  of  r  have  been 
measured,  the  thermal  conductivity  may  be  calc\ilated  by  integrating 
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Eq.  (20)  "between  the  wall  and  the  center  of  the  arc,  giving 


The  temperature  gradient  is  measured  spectroscopically  along  a 
line  perpendicular  to  the  axis  of  the  arc.  Most  authors  have  used  inten¬ 
sity  measurement  of  line  radiation  assuming  the.  transition  probability  is 
known  and  checked  the  temperature  by  measuring  the  continuum  intensity. 
The  temperature  measurement  relies  on  the  gas  being  optically  thin  at 
the  wavelength  used  for  the  temperature  determination.  The  population 
of  electronic  states  and  the  ionization  must  also  be  in  equilibrium  with 
the  translational  temperature  of  the  heavy  particles.  The  rate  at  which 
the  electronic  states  are  populated  by  electron- heavy  particle  collisions 
must  be  much  larger  than  the  rate  at  which  depopulation  by  spontaneous 
radiation  occurs  in  order  to  reconcile  the  optically  thin  criterion  with 
the  equilibrixim  population  of  electronic  states.  The  first  problem  is  to 
obtain  the  intensity  of  the  radiation  as  a  function  of  r  from  the  integrated 
intensity  measurements.  Solution  of  this  problem  depends  on  the  cylin¬ 
drical  symmetry  of  the  arc  and  is  discussed  below  in  the  section  on 
radiation.. 

The  spectroscope  sees  the  radiation  from  the  volume  element 
along  a  cord  of  the  arc  as  shown  in  Figure  26.  The  optics  are  held 
stationary  and  the  arc  is  moved  perpendicular  to  the  arc  axis  and  another 
measurement  of  radiative  intensity  is  made.  The  intensity  from  such 
volume  elements  as  a  function  of  distance  from  the  center  of  the  arc 
S(x)  is  thus  obtained.  The  desired  quantity  is  the  radiative  intensity 
from  a  volume  element  around  a  point  at  a  radius  r  from  the  center  of 
the  arc,  I(r).  An  integral  relation  between  the  two  quantities  may  be 
derived  by  considering  the  volume  element  along  the  cord  to  be  made 
up  of  cubic  elements  at  different  radii.  From  the  geometry  of  Figure 
26  we  have 

S(x)  =  J  I(r)  dy  .  (22) 

Relating  the  distance  along  the  cord  y  to  x  and  r  an  integral  equation 
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relating  S(x)  to  l(r)  is  obtained 


S(x)  = 


7^ 


(23) 


Since  S(x)  and  x(r)  are  known  functions,  I{r)  may  be  determined  analyt¬ 
ically  using  Abel' s  transformation  when  radial  symmetry  is  assumed. 
The  inversion  process  is  discussed  at  length  by  several  authors  (Refs. 

50,  80)  and  leads  typically  to  about  a  2%  error  in  the  temperature  profile. 
The  most  difficult  experimental  problem  is  maintaining  the  absolute 
calibration  of  the  intensity  measurement.  This  can  be  done  by,  for  in¬ 
stance,  comparison  with  a  standard  source  such  as  a  calibrated  carbon 
filament. 


The  second  problem  is  to  relate  the  intensity  radiated  by  a  volume 
element  to  its  temperature.  The  type  of  radiation  (line,  molecular  band, 
continuum)  used  for  the  temperature  measurement  depends  largely  on  the 
temperature  and  species  present.  The  radiative  intensity  from  line 
radiation  is  better  known  than  other  forms  of  radiation.  Radiative  in¬ 
tensity  versus  temperature  may  also  be  constructed  for  other  types  of 
radiation  such  as  continuum  and  molecular  radiation.  However,  theory 
may  be  relied  on  only  for  the  temperature  dependence  of  other  types  of 
radiation,  and  the  absolute  value  shovild  be  checked  by  measuring  atomic 
or  ionic  lines  in  the  same  frequency  range. 


The  absolute  intensity  per  unit  volume  and  solid  angle  of  an 
atomic  line  is  given  by  (Ref.  80) 
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where  C  is  the  speed  of  light,  f  the  length  of  the  source  (volxome  ele¬ 
ment),  h  the  Planck  constant,  the  number  density  of  atoms,  Q  the 
partition  function  of  the  atomic  specie,  the  statistical  weight  of  the 
upper  state  of  the  transition,  the  energy  of  the  upper  level,  V  the 
wave  niomber  of  the  line,  and  A^^^  the  transition  probability. 

The  niunber  of  atoms  and  the  partition  function  are  known  rather 
well  (~2%)  for  an  equilibrium  gas  from  the  usual  thermodynamic 
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calculations.  The  transition  probabilities  have  large  uncertainties. 
However,  large  uncertainties  in  the  transition  probabilities  {~20%)  lead 
to  small  uncertainties  in  the  gas  temperature  (Ref.  50).  Equation  (24) 
may  be  used  to  obtain  a  table  of  intensity  per  unit  volume  versus  tem¬ 
perature. 

The  radiative  loss  term  in  Eq.  (20)  is  usually  determined  by  in¬ 
tegrating  the  measured  radiation  from  the  arc,  the  assumption  again 
being  that  the  gas  is  optically  thin.  Early  arc  measurements  indicated 
that  over  much  of  the  temperature  range  (N21  4000-7000°K;  Ar,  6OOO- 
10,  000  K)  the  radiation  is  not  large  and  the  optically  thin  assumption 
is  adequate.  However,  recently  the  group  at  Avco  showed  that  apply¬ 
ing  appropriate  radiation  corrections  to  Maecker*  s  data  removed  the 
anomalous  current  dependence  of  the  nitrogen  thermal  conductivity. 

In  addition,  the  experimental  thermal  conductivities  from  Maecker*  s 
experiment  are  in  surprisingly  good  agreement  with  recent  transport 
property  calculations.  The  large  scatter  in  the  Maecker  arc  runs  occurs 
near  the  dissociation  peak  (~7000  K).  Both  the  Maecker  and  Avco  data 
also  show  a  marked  current  dependence  (~  a  factor  of  three  between  250 
amp  runs  and  17  amp  runs)  above  8000  K  as  a  free-free  and  free-bound 
continuum  becomes  important.  Another  approach  was  taken  by  Asinovskii 
and  Kirillin  (Ref.  49)  where  the  thermal  conductivity  of  argon  was  deter¬ 
mined  in  arcs  of  different  diameters.  As  the  arc  diameter  was  reduced 
the  thermal  conductivity  fell  closer  to  the  theoretical  value.  The  thermal 
conductivity  versus  arc  radius  from  this  experiment  is  shown  in  Figure 
27  along  with  the  theoretical  value  (see  argon  data  review.  Figure  7). 

Also  shown  in  Figure  27  are  the  argon  data  of  Avco  and  that  of  Emmons 
and  Knopp.  The  high  current  data  of  the  first  two  investigators  lies  much 
higher  than  the  low  current  data  of  Emmons  and  Knopp,  indicating  the  im¬ 
portance  of  radiative  transport. 

High  arc  currents  mean  high  centerline  temperatures.  The  vac¬ 
uum  ultraviolet  becomes  more  important  at  the  high  temperatures  because 
the  continuum  cutoff  is  higher  and  the  higher  energy  levels  radiate  more. 
The  radiative  transport  term  Q  in  Eq.  (20)  then  includes  absorption. 

Since  the  temperature  is  known  as  a  function  of  position,  the  radiative 
intensity  may  be  calculated  at  each  position  in  the  arc  using  the  equation 
of  transfer 
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where  r'  is  the  radial  distance  from  the  element  of  fluid  under 


26 


AEDC-TR.68>105 


consideration  to  the  source  of  the  radiation,  is  the  absorption  coef¬ 
ficient,  is  the  Planck  function  and  is  the  optical  path  between  the 
point  at  which  I(r)  is  being  calculated  and  the  source  point  r*  , 


(26) 


The  absorbed  radiation  in  the  voliime  element  is  the  previously  calcu¬ 
lated  intensity  (Eq.  25)  times  the  spectral  absorption  coefficient  inte¬ 
grated  over  all  solid  angles  and  all  frequencies.  The  net  flux  of  energy 
Q  which  must  be  inserted  into  the  energy  equation  (20)  is  the  difference 
between  the  absorbed  energy  and  the  radiated  energy 


Q{r)  = 


I  d^  -  4  a 

V  V 


(27) 


When  the  usual  opticaUy  thin  approximation  is  applicable  a  is  small 
and  recalling  that 


I  =  a  B  , 

V  V  V 


(28) 


Q(r)  in  Eq.  (20)  reduces  to  the  usual  form, 

Q(r)  =  Att  r  dv  (29) 

since  the  first  term  is  second  order  in  a^. 

The  major  problem  is  to  obtain  values  of  spectral  absorption 
coefficient  and  optical  path  at  the  far  ultraviolet  wavelengths. 
Theoretical  calculations  are  available  from  many  sources.  Experi¬ 
mentally,  however,  plasma  radiation  has  not  been  measured  below 
2000A.  The  major  reason  is  that  the  vacuum  ultraviolet  is  absorbed 
in  the  cool  outer  layer  of  the  arc.  Indications  are  that  continuum  ultra¬ 
violet  should  be  significant  down  to  about  800A  (Refs.  81,  82,  83,  84,  85,  86) 
in  the  gases  and  temperatures  of  current  interest.  Calculations  of 
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radiative  losses  in  acoustic  waves  indicate  that  ultrasonic  absorption 
may  provide  information  on  the  vacuum  ultraviolet  (Refs.  9,  87,  88,  89). 
The  acoustic  technique  has  the  advantage  of  reaching  into  the  hot  gas 
thus  avoiding  the  boundary  layer  problems  of  other  techniques. 

Finally,  the  electrical  heating  term  aE  in  Eq.  (20)  has  to  be 
evaluated.  The  electrical  field  strength  may  be  measured  by  two  tech¬ 
niques.  First,  the  length  of  the  arc  is  increased  by  moving  the  elec¬ 
trodes  a  known  distance,  keeping  the  total  arc  current  constant.  The 
slope  of  the  resulting  voltage  drop  versus  arc  length  is  the  electrical 
intensity.  Accurate  measurements  require  that  the  voltage  drop  in  the 
region  of  cathode  and  anode  be  a  reasonably  small  faction  of  the  voltage 
drop  in  the  central  portions  of  the  arc.  A  second  method  is  used  in  arcs 
made  up  of  copper  plates  alternated  with  plates  of  insulating  material. 
The  voltage  of  successive  plates  near  the  center  of  the  arc  are  plotted 
versus  position  and  the  slope  provides  the  field  intensity.  Typical 
errors  reported  for  field  intensity  measurements  are  6%. 

The  electrical  conductivity  may  in  principal  be  determined  from 
the  known  temperature  distribution,  field  intensity,  and  the  measured 
total  current.  The  electrical  conductivity  as  a  fxinction  of  temperature 
(Ref.  49)  may  be  obtained  by  inverting  the  integral  equation 

R 

ar  dr  (30) 


where  R  is  the  radius  of  the  arc  and  O’  and  r  are  known  functions  of  T. 

In  practice  nonequilibrium  and  contact  potentials  in  the  cooler  regions 
of  the  arc  may  lead  to  errors  and  calculated  values  are  preferable.  The 
most  recent  calculation  and  probably  the  most  reliable  are  given  by 
Devoto  (Ref.  21)  for  argon.  Argon  electrical  conductivities  have  been 
reviewed  by  others.  For  nitrogen  only  theoretical  electrical  conduc¬ 
tivities  are  available.  Electrical  conductivity  is  of  interest  in  its  own 
right  as  a  transport  property. 

The  constricted  arc  experiment  has  produced  some  of  the  best 
thermal  conductivity  measurements  obtained  to  date  at  extremely  high 
temperatures.  Radiative  transport  in  the  far  ultraviolet  is  important 
at  about  50%  dissociation  and  once  ionization  becomes  significant. 
Knowledge  of  the  absorption  coefficient  is  partictilarly  necessary  at 
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temperatures  above  1300°K  and  for  work  at  higher  pressures.  In  the 
high  temperature  (13,  000°K  <  T  for  argon)  and  pressure  regions 
(1  atm  <  p)  the  radiative  term  should  be  evaluated  using  the  optically 
thick  equation  of  transfer.  Theoretical  electrical  conductivities  seem 
adequate  (within  ~6%)  for  argon.  However,  independent  electrical  con¬ 
ductivity  measurements  should  be  obtained  in  nitrogen  and  other  gases 
before  the  constricted  arc  is  used  to  measure  thermal  conductivity. 

The  effect  of  pressure  on  electrical  conductivity  must  also  be  investi¬ 
gated  before  the  constricted  arc  is  extended  to  higher  pressures. 

3.3  HEAT  TRANSFER 

The  heat  transfer  to  the  end  wall  after  the  reflection  of  a  normal 
shock  in  the  shock  tube  has  been  suggested  as  a  means  of  obtaining  ther¬ 
mal  conductivity  at  high  temperatures  (Refs.  90,  91).  The  reflection  of 
the  incident  shock  in  the  first  approximation  creates  a  uniform  sample  of 
high  temperature  gas  extending  about  a  foot  between  the  reflected  shock 
front  and  the  end  wall.  In  the  vicinity  of  the  end  wall  a  thermal  boundary 
layer  is  set  up.  Since  the  boundary  layer  is  thin  the  pressure  is  con¬ 
stant  throughout  and  equal  to  the  free  stream  pressure.  The  gas  tem¬ 
perature  and  density  in  the  end  wall  boundary  layer  varies  between  the 
free  stream  temperature  of  the  reflected  shock  region  and  the  wall  tem¬ 
perature.  The  end  wall  geometry  is  shown  in  Figure  28.  The  isotherms 
(which  are  also  the  isopycnals)  are  shown  as  lines  parallel  to  the  end 
wall  and  the  reflected  shock  front.  The  close  spacing  of  the  isotherms 
near  the  wall  indicates  a  steepening  of  the  temperature  gradient  near  the 
wall.  Because  of  the  large  mass  of  the  shock  tube  the  end  wall  experi¬ 
ences  only  a  slight  temperature  rise  as  it  cools  the  boundary  layer  gases. 

The  experimental  measurement  of  heat  flux  to  the  end  wall  of  the 
shock  tube  is  very  simple.  The  heat  transfer  gauges  are  1  cm  thick 
Pyrex  or  quartz  blanks  with  a  thin  film  platinum  thermometer  deposited 
on  the  surface  as  shown  in  Figure  28.  The  surface  temperature  of  the 
Pyrex  blank  is  an  integral  function  of  the  heat  flux  versus  time  to  the 
surface  of  the  blank.  Therefore,  the  thin  film  thermometer  response 
is  an  integral  function  of  the  heat  flux  versus  time.  In  the  usual  mode 
of  operation  the  thin  film  thermometer  is  made  sufficiently  thin  so  that 
it  measures  the  Pyrex  surface  temperature  without  appreciably  perturb¬ 
ing  the  temperature  versus  time  history  in  the  Pyrex  slab. 

The  ideal  end  wall  response  is  a  step  increase  in  temperature 
(voltage)upon  reflection  of  the  shock.  Boundary  layer  theory  shows  that 
the  heat  flvix  to  the  end  wall  of  a  shock  tube  should  decrease  with  the 
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square  root  of  time  as  the  boundary  layer  grows  (Refs.  92,  93).  Heat 
transfer  theory  (Ref.  94)  indicates  that  the  response  of  an  infinite  slat 
of  material  subjected  to  such  a  heat  flux  is  a  step  rise  in  surface  tem¬ 
perature.  The  temperature  rise  is  given  by 
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This  in  turn  is  related  to  the  enthalpy  gradient  at  the  wall  by  Eq,  (32). 
The  quantities  p  Cp5v,  are  the  properties  of  the  glass  blank.  In  the  first 
approximation  the  thermal  properties  of  the  Pyrex  are  constant  since 
the  temperature  rise  at  the  wall  is  usually  of  the  order  of  10  K.  A  cor¬ 
rection  for  the  temperature  dependence  of  p  c^A  has  been  worked  out  by 
Hartunian  and  Varwig  (Ref.  95). 

The  heat  flux  from  the  end  wall  boundary  layer  of  the  shock  tube 
is  susceptible  to  simple  analysis.  The  flow  in  the  boundary  layer  is 
described  (after  a  small  transient  effect)  by  the  one  dimensional  Navier- 
Stokes  conservation  equations  of  mass,  energy  momentum,  and  by  the 
equation  of  state.  The  equations  of  motion  are  further  simplified 
because  the  pressure  is  constant  in  the  boundary  layer  and  viscous 
dissipation  may  be  neglected.  However,  convection  terms  must  be  re¬ 
tained  because  of  thermal  expansion.  The  early  attempts  (Refs.  96,  97) 
to  measure  thermal  conductivity  of  polyatomic  gases  on  the  end  wall  of 
the  shock  tube  were  in  error  since  they  neglected  convection  and  disso¬ 
ciation  relaxation. 


One  of  the  most  frequently  misunderstood  concepts  is  the  depend¬ 
ence  of  the  flow  at  the  wall  on  the  conditions  throughout  the  boundary 
layer.  The  heat  flux  to  the  wall  may  be  written  in  the  following  way; 


Q  =  -\ 
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where  Q  is  the  heat  flux  and  k  ^  is  the  thermal  conductivity  at  the  wall 
temperature.  The  temperature  gradient  at  the  wall  depends  on  the 
solution  to  the  boundary  layer  equations  subject  to  the  free  stream  and 
wall  boundary  conditions.  The  temperature  gradient  at  the  wall  there¬ 
fore  reflects  the  behavior  of  thermal  conductivity  throughout  the  boundary 
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layer. 


The  energy  equation  (8)  for  the  end  wall  boundary  layer  may  be 
written  as 


where  the  distance  from  the  end  wall  y  has  been  replaced  by 


1  =  7^  f  P  dy  (34) 

as  the  independent  variable  to  reduce  the  energy  equation  to  the  incom¬ 
pressible  form,  p  is  the  density,  Cp  the  specific  heat  in  the  boundary 
layer  and  \  the  thermal  conductivity.  The  enthalpy 


w 


(35) 


replaces  temperature  as  the  dependent  variable.  The  quantity  is 

enthalpy  dependent  coefficient  of  Eq.  (33).  P 

Once  is  a  known  function  of  h  Eq.  (33)  is  a  second  order  dif- 

ferential  equation  which  has  a  unique  solution 


h  =  h  (n)  =  h  (x,t) 


(36) 


if  two  boundary  conditions  are  specified.  The  boundary  conditions  are 
that  the  enthalpy  at  the  wall  is  zero  and  the  enthalpy  at  the  free  stream 
(in  the  reflected  shock  region)  is  known  at  reflected  shock  pressures 
and  temperatures. 

For  any  gas  where  is  known  as  a  function  of  the  enthalpy  and 
pressure  (which  are  equivalent  to  temperature  and  pressure),  Eq.  (33) 
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has  a  unique  solution,  h  (x,  t),  which  may  be  differentiated  to  calculate 
the  heat  flux  to  the  wall  according  to  Eq.  (32).  This  process  may  be 
turned  around  to  obtain  the  temperature  dependence  of  £A  from  the 

measurements  of  heat  flux  and  reflected  shock  conditions  in  the  shock 
tube.  Because  the  heat  flux  depends  on  the  functional  relationship  be¬ 
tween  thermal  conditivity  and  temperature  rather  than  on  the  value  of 
the  thermal  conductivity  at  one  temperature,  we  say  that  the  heat  flux 
is  a  functional  of  thermal  conductivity. 

An  integral  relation  to  replace  the  differential  equation  (33)  which 
was  derived  by  Jepson  (Ref.  92)  is  discussed  in  Appendix  I.  This  rela¬ 
tion  provides  the  heat  flux  in  terms  of  the  transport  properties  and  the 
enthalpy  range  of  the  boundary  layer.  The  expression  is 

L  ^  (”) 

w  w  P 


Equation  (37)  has  been  investigated  as  a  technique  for  determining  the 
temperature  dependence  of  P^L  .  It  was  found  that  extremely  accurate 

^P 


(Ref.  98)  heat  transfer  measurements  over  a  wide  range  of  temperatures 
were  necessary  to  differentiate  between  various  possible  functional  de¬ 
pendencies  for  which  would  give  similar  values  for  Q. 


However,  both  theoretically  and  experimentally,  the  end  wall  ex¬ 
periment  is  suited  to  determine  the  density  dependence  of  thermal  con¬ 
ductivity.  The  theory  takes  advantage  of  the  fact  that  the  excess  thermal 
conductivity  is  a  function  of  density  alone.  In  Sengers'  notation  (Ref.  68) 
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where  k  ^  is  the  low  pressure  thermal  conductivity,  p  is  the  density  and 
Aj^  ,  are  the  virial  coefficients  of  the  thermal  conductivity.  Since  the 
virial  coefficients  are  essentially  independent  of  temperature,  the  heat 
flux  then  becomes 
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The  first  term  is  the  heat  flux  which  may  be  calculated  from  the  pre¬ 
sumably  known  values  of  the  low  pressure  thermal  conductivity,  \  q, 
density,  p  ,  and  specific  heat,  Cp.  The  other  integral  is  a  known 
thermodynamic  quantity  except  for  multiplicative  constants  indicated 
by  Equations  (37)  and  (38). 

Writing  the  above  equation  explicitly  in  terms  of  the  derived  quan¬ 
tities  Ay^  and  from  Eq.  (38) 

=  Q  ^  +  Q,^  A.  +  B.  (40) 

O  X  \  u  K 


where  Qj  and  Q2  are  integrals  over  functions  of  density.  Measurements 
of  Q  at  at  least  two  temperatures  will  allow  Aj^  and  Bj^  to  be  calculated. 
Ay^  and  Bj^  may  not  correspond  to  the  theoretical  virial  coefficients  due 
to  the  fact  that  it  is  difficult  to  separate  the  effect  of  higher  order  terms 
in  the  density  expansion.  However,  the  density  dependence  of  \  -\q  may 
be  established  to  about  1%. 

For  the  pressure  dependence  of  thermal- conductivity  we  are  inter¬ 
ested  in  free  stream  temperatures  of  less  than  2000  K.  Ionization  and 
dissociation  relaxation  do  not  occur  at  these  temperatures.  The  effects 
of  shock  attenuation  (Ref.  8)  may  be  avoided  by  making  measurements  a 
few  psecs  after  the  reflection  of  the  incident  shock.  The  thermodynamic 
state  is  determined  by  wave  speed  measurements  close  to  the  end  wall. 
This  should  be  corroborated  with  a  pressure  measurement  using  a  fast 
response  pressure  gauge. 

The  calibration  of  the  heat  transfer  gauge  must  be  carefully  mon¬ 
itored  for  accurate  results.  The  end  wall  thin  film  gauges  tend  to  pit  and 
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wear  (Ref.  8).  A  measurement  of  gauge  resistance  after  each  shock 
is  necessary  at  the  high  Mach  number  and  presumably  at  high  pres¬ 
sures.  Determination  of  the  thermal  coefficient  of  resistivity  should 
be  carried  out  for  each  gauge.  The  gauge-  characteristics  have  to  be 
calibrated  after  every  four  or  five  shocks.  The  easiest  method  of  cal¬ 
ibration  is  firing  low  Mach  number  shocks  in  a  standard  gas  such  as 
argon  or  nitrogen.  At  the  extremes  of  pressure,  p  ~  10'*  psi,  the  gauge 
may  be  destroyed  during  the  measurement. 

The  heat  transfer  technique  provides  a  simple  method  to  obtain 
information  concerning  thermal  conductivity  at  high  temperature  and 
pressure.  The  measured  quantity,  heat  transfer  to  the  wall,  depends 
on  the  thermal  conductivity  throughout  the  boundary  layer.  Therefore, 
extremely  accurate  heat  transfer  measurements  (~0.  1%)  are  necessary. 
In  addition,  care  must  be  taken  to  avoid  relaxation  effects  such  as  ioniza¬ 
tion  and  dissociation.  At  low  densities  the  heat  transfer  method  mea¬ 
sures  the  temperature  dependence  of  thermal  conductivity.  However, 
the  technique  is  very  insensitive  to  the  thermal  conductivity  at  the  higher 
temperatures  in  the  boundary  layer. 

The  heat  transfer  technique  is  best  suited  to  measure  variations 
in  thermal  conductivity  which  are  systematic  with  temperature.  For 
instance,  the  heat  transfer  technique  can  choose  the  best  parameters 
for  a  given  potential  function.  Another  area  where  heat  transfer  mea¬ 
surements  may  contribute  is  in  establishing  the  mixing  rules  for  the 
thermal  conductivity  at  high  temperatures  and  pressures.  Finally,  as 
discussed  above,  the  heat  transfer  technique  is  particularly  suitable  for 
the  determination  of  the  density  dependence  of  thermal  conductivity. 

3.4  SCHLIEREN  TECHNIQUE 

Optical  techniques  may  be  used  to  measure  density  and  density 
gradients  in  the  end  wall  boundary  layer  in  the  shock  tube  (Ref.  8). 

This  information  together  with  the  energy  equation  for  the  end  wall 
boundary  layer  may  be  used  to  obtain  thermal  conductivity.  Measure¬ 
ments  have  been  carried  out  by  Smeets  (Ref.  99)  in  air  up  to  6000  K. 
During  the  past  year  the  schlieren  system  shown  in  Figure  29  was  set 
up  on  an  end  wall  of  the  shock  tube  for  thermal  conductivity  measure¬ 
ments. 

A  schematic  of  the  system  is  shown  in  Figure  30.  A  line  source 
of  light  is  made  by  mounting  a  mercury  flash  tube  A  behind  a  source 
slit  S]^,  located  at  the  focus  of  the  collimating  lens  B.  Parallel  light 


34 


AEDC.TR.68.105 


then  passes  through  the  reflected  shock  region.  The  arrow  (Figure  30) 
pointing  away  from  the  end  wall  denotes  the  object  which  is  imaged  at  the 
film  plane.  The  boundary  payer  extends  from  the  foot  of  the  arrow  to  the 
head  where  the  free  stream  begins. 

Rays  which  travel  through  the  head  of  the  arrow  remain  parallel 
because  the  index  of  refraction  is  constant  in  the  free  stream  and  cor¬ 
responds  to  Pg.  Rays  which  pass  closer  to  the  end  wall  travel  slower 
than  rays  passing  through  the  free  stream  because  the  density  near  the 
wall  is  greater  than  the  free  stream  density.  Therefore  the  wave  front 
D  passing  through  the  region  close  to  the  end  wall  is  bent  toward  the  wall. 

The  parallel  light  from  the  free  stream  region  is  focused  by  the 
mirror  C  at  the  focal  point  of  C.  A  knife  edge  S2  is  put  into  the  field 
so  as  to  block  the  undeviated  light.  The  free  stream  region  then  appears 
black  on  the  photographic  plate.  On  the  other  hand  the  light  traveling 
through  the  boundary  layer  is  bent  so  that  it  is  no  longer  parallel  to  the 
undisturbed  light.  Consequently,  the  rays  from  the  image  closer  to  the 
wall  are  focused  at  a  point  above  the  knife  edge  S^.  The  boundary  layer 
region  therefore  appears  light  on  the  photograph.  This  is  illustrated  in 
the  schlieren  photograph  of  the  boundary  layer  taken  (Figure  31)  400  psec 
behind  a  reflected  shock  in  argon.  A  densitometer  trace  of  the  contrast 
on  the  negative  of  the  photograph  ideally  will  provide  the  density  gradient 
in  the  boundary  layer.  Putting  the  density  gradient  in  the  linearized 
energy  equation  (Ref,  8)  for  the  end  wall  boundary  layer  equations  allows 
the  thermal  conductivity  to  be  determined 


C  ^  ^ 

^^2F  p 


(41) 


where  C  is  the  measured  contrast  on  the  schlieren  photograph. 

Quantitative  measurements  with  the  usual  schlieren  system  has 
proved  extremely  difficult  because  of  effects  due  to  diffraction,  light 
reflection  from  the  end  wall,  and  uneven  illumination.  These  difficulties 
are  largely  eliminated  by  using  schlieren  interferometry  (Refs.  100,  101, 
102).  The  fringes  are  produced  by  replacing  the  end  slit  with  a  wire 
of  small  diameter  D  subject  to  the  condition 


D  <  1.  2  \  /tt  d 


(42) 
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where  f  is  the  focal  length  of  the  mirror,  X  is  the  wavelength  of  the 
source,  and  d  is  the  thickness  of  the  boundary  layer.  The  schlieren 
interferometer  described  above  responds  to  density,  not  density  gradi¬ 
ents,  and  two  differentiations  are  necessary  to  determine  thermal  con¬ 
ductivity  from  the  end  wall  boundary  layer  equation. 

An  interferometric  method  in  which  the  fringe  shift  is  proportional 
to  the  density  gradient  has  been  used  by  Smeets  to  determine  the  thermal 
conductivity  of  air  from  measurements  of  the  end  wall  boundary  layer. 

The  experimental  arrangement  is  similar  to  that  pictures  in  Figure  30. 
The  essential  new  features  are  the  addition  of  polarizers  and  Wollaston 
prisms  between  Sj  and  the  collimating  lens  B  as  well  as  between  the 
focal  point  {S2)  of  the  schlieren  lens  C,  and  the  camera  lens  G.  The 
Wollaston  prisms  are  placed  in  symmetrical  positions  in  the  source 
and  image  parts  of  the  optical  system.  Further  details  which  improve 
the  performance  of  the  optics  are  given  in  the  original  paper  (Ref.  103). 
Smeets  reduced  his  interferometer  measurements  by  calculating  the  ex¬ 
pected  fringe  shift  from  a  linearized  approximation  to  Eq.  (33).  The 
thermal  conductivity  at  the  free  stream  was  determined  by  making  the 
expected  fringe  shift  agree  with  the  measured  fringe  shift.  Data  reduc¬ 
tion  for  this  type  of  experiment  is  discussed  further  in  Appendix  I, 

Optical  measurements  on  the  end  wall  of  the  shock  tube  are  one  of 
the  few  ways  of  obtaining  high  temperature  thermal  conductivities.  In 
addition,  the  pressure  in  the  reflected  region  of  the  shock  tube  may  be 
varied  easily.  However,  measurements  must  be  made  at  long  times 
(T  >  400  |JLsec)  and  low  pressures  so  that  the  boundary  layer  is  suffi- 
ently  thick  to  be  resolved  optically.  Long  times  are  also  necessary  to 
a,llow  relaxation  due  to  dissociation  and  ionization  to  occur  and  the  bound¬ 
ary  layer  to  return  to  equilibrium.  The  pressure  and  temperature  in  the 
reflected  region  can  vary  up  to  20%  in  times  of  the  order  of  100  [xsec  due 
to  attenuation  and  boundary  layer  interaction.  Therefore  the  pressure 
and  temperature  should  be  measured  simultaneously  with  the  optical 
measurements. 


3.5  ULTRASONIC  TECHNIQUE 

The  ultrasonic  technique  has  recently  been  developed  for  the  de¬ 
termination  of  transport  properties  of  high  temperature  gases.  The 
absorption  of  ultrasonic  waves  in  a  gas  is  well  understood  and  depends 
on  many  factors.  In  monatomic  gases  at  low  temperatures  the  sound 
absorption  depends  on  thermal  conductivity  and  viscosity  and  is  given  by 
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where  u  is  the  angular  frequency  of  the  sound  wave  and  the  other  quan¬ 
tities  are  previously  defined.  Since  viscosity  is  much  easier  to  measure 
than  thermal  conductivity,  measured  values  of  viscosity  and  sound  ab¬ 
sorption  will  provide  a  useful  technique  for  the  measurement  of  thermal 
conductivity.  This  technique  has  been  applied  in  helium  and  argon  up  to 
temperatures  of  about  11,  000  K. 

Sound  absorption  in  mixtures  of  monatomic  gases  (Refs.  9,  104) 
depends  on  the  diffusion  coefficients  as  well  as  thermal  conductivity  and 
viscosity  of  the  mixture.  The  absorption  is  then  given  by 
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These  new  quantities  are  the  thermal  diffusion  coefficient  a  ij. ,  the  mole 
fractions  xj^,  x^,  and  the  atomic  masses  m^^  and  m2  (mj  >  m2).  The 
size  of  the  compared  to  depends  largely  on  the  size  of  the  term 
in  brackets  and  the  mole  fractions.  When  the  ajj  is  sufficiently  large 
the  ultrasonic  technique  is  probably  the  best  technique  for  measuring 
diffusion  coefficients.  Ultrasonic  measurement  of  diffusion  coefficients 
have  been  made  by  several  authors  (Refs.  105,  106,28),  the  highest 
temperature  being  5000  K  in  the  reflected  region  of  the  shock  tube. 

Polyatomic  gases  and  inert  gases  at  high  temperatures  are  more 
complicated  because  of  relaxation  processes  such  as  vibration,  rotation, 
electronic  excitation  and  chemical  reactions  (Ref.  9).  At  sufficiently 
low  temperatures  (T  <  2000°K  for  O2,  T  <  5000  K  for  N2)  there  are  no 
chemical  degrees  of  freedom.  Vibration  and, rotation  are  the  most  com¬ 
mon  internal  degrees  of  freedom.  The  absorption  due  to  a  single  inter¬ 
nal  degree  of  freedom  is  given  by  (Ref.  107) 
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where  Cp  and  are  the  total  or  equilibrium  of  specific  heats  and  is 
the  specific  heat  of  the  internal  degree  of  freedom.  The  relaxation  time 
T*  is  related  to  the  true  relaxation  time  of  the  internal  mode  T  by  the 
relation 


C  -C. 
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(47) 


and  the  true  relaxation  time  is  defined  by  the  equation  of  change  (Ref.  15) 
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The  internal  modes  also  cause  a  velocity  dispersion.  The  ratio  of 
the  low  frequency  limit  sound  speed  (c  )  to  the  sound  speed  at  the  operat¬ 
ing  frequency  is 
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where  the  low  frequency  sound  speed  is 
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The  physical  content  of  these  equations  is  that  the  energy  in  the 
internal  modes  do  not  follow  the  energy  in  the  translational  modes  for 
all  frequencies.  The  sound  wave  consists  of  periodic  fluctuations  of 
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the  kinetic  energy  of  the  gas  particles.  In  the  process  the  internal  mode 
absorbs  energy  from  translation  in  the  high  temperature  regions  and 
releases  energy  to  translation  in  the  low  temperature  regions.  The  rate 
of  the  energy  exchange  between  an  internal  mode  with  translation  is 
smaller  than  the  rate  of  translation- translation  exchange.  The  exchange 
process  may  be  described  by  the  time  required  for  the  internal  modes  to 
come  isothermally  (Ref.  107)  to  equilibrium  with  translation  after  a  step 
in  translational  temperature.  This  is  identical  to  the  relaxation  time  de¬ 
fined  by  Eq.  (48). 

When  the  period  of  the  wave  is  much  greater  than  the  relaxation 
time  for  energy  exchange,  that  is,  at  low  frequencies  uT  <<  1,  the 
sound  speed  and  absorption  contain  contributions  from  the  internal 
modes.  As  the  frequency  of  the  sound  wave  is  increased  the  energy  in 
the  internal  modes  begins  to  fall  off  because  it  becomes  more  difficult 
for  them  to  adjust  to  the  changes  in  translational  temperature.  Finally, 
at  frequencies  so  high  (wT  >>  1)  that  the  period  of  the  sound  wave  is 
much  shorter  than  the  relaxation  time  the  internal  mode  exchanges  very 
little  energy  with  the  sound  wave.  The  sound  wave  then  propagates  as 
though  the  slow  internal  mode  were  not  there  and  the  internal  mode  is 
said  to  be  frozen  out. 

The  sound  absorption  in  a  gas  with  internal  modes  is  made  up  of 
the  classical  contributions  given  by  Equations  (43)  and  (45)  with  the  ab¬ 
sorption  due  to  relaxation  as  given  by  (46)  added  on.  Although  relaxation 
times  are  of  interest  in  themselves  the  primary  interest  is  in  transport 
properties.  Therefore,  the  absorption  due  to  internal  modes  must  be 
subtracted  from  the  total  absorption.'  There  are  several  methods  for 
handling  these  problems.  The  simplest  technique  is  simply  to  operate 
at  sufficiently  high  frequencies  that  the  internal  mode  is  frozen  out. 

This  may  be  done  for  vibration,  dissociation  and  ionization  with  ultra¬ 
sonic  waves  in  the  megahertz  region  up  to  10,  000  K.  Alternatively  the 
pressure  may  be  varied. 

The  second  approach  is  to  use  auxiliary  experiments  to  define  the 
relaxation  times  and  calculate  the  absorption  due  to  internal  modes  using 
Eq.  (46).  Rotation  and  electronic  excitation  must  be  handled  in  this  fash¬ 
ion.  For  internal  modes  which  have  relaxation  times  longer  than  roughly 
10“ ®  sec,  ultrasonic  absorption  vs  frequency  will  provide  the  necessary 
information.  There  are  many  other  techniques  such'  as  thermal  trans¬ 
piration  (Ref.  30),  electron  beam  fluorescence  (Ref.  29),  as  well  as 
various  shock  tube  experiments  which  can  provide  relaxation  times  where 
necessary. 
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The  experimental  problems  associated  with  the  ultrasonic  tech¬ 
nique  depend  upon  the  pres  sure- temperature  range  of  the  measurements. 
At  temperatures  which  may  be  obtained  in  conventional  furnaces  (2000 
to  2500  K),  the  transducers  have  to  be  protected  by  buffer  rods  of  high 
temperature  materials  or  water- cooling  must  be  employed.  Recently 
developed  piezoelectric  materials  can  be  used  at  temperatures  up  to 
1000  C  (Ref.  108).  At  least  two  acoustic  paths  must  be  used  to  cancel 
Out  the  effects  of  thermal  boundary  layers.  At  higher  temperatures 
momentary  contact  must  be  employed. 

The  main  problem  at  high  pressures  is  the  necessity  of  using  high 
frequency  sound  waves  since  the  absorption  in  general  decreases  with  an 
increase  in  pressure.  Generally,  a  change  of  about  10  db  is  necessary 
to  obtain  a  good  absorption  measurement.  Typically,  1-3  MHz  give  good 
results  at  one  atmosphere  and  temperatures  of  the  order  of  10,  000  K. 
Since  absorption  increases  with  ,  increasing  the  frequency  by  a  factor 
^  P 

of  vp  maintains  a  constant  sound  absorption  as  the  pressure  is  increased. 
Therefore,  extension  of  the  ultrasonic  technique  to  lO'^  atmosphere  re¬ 
quires  frequencies  of  the  order  of  100  MHz. 

Arc  heated  gases  have  been  probed  ultrasonically  between  5000°K 
and  20,  000  K.  The  lower  temperature  limit  is  imposed  because  a  cer¬ 
tain  amount  of  ionization  is  necessary  for  arc  heating  (Ref.  9).  It  is 
possible  to  probe  the  exhaust  of  arc  heated  gases  sufficiently  down  stream 
to  obtain  gas  temperatures  below  5000  K.  However,  arc  exhausts  tend  to 
be  very  turbulent  leading  to  scattering  of  the  ultrasonic  wave  and  acoustic 
noise  from  pressure  fluctuations.  High  pressure  arcs  have  high  radiation 
losses  which  make  the  power  supply  and  cooling  problem  more  difficult. 

In  addition,  at  high  pressures  radiation  becomes  a  significant  transport 
mechanism. 

The  most  promising  gas  source  in  the  intermediate  temperature 
range  1000  K  to  10,  000  K  is  the  reflected  region  of  the  shock  tube. 
Several' gases  have  been  probed  in  the  shock  tube  and  agreement  with 
ultrasonic  measurements  in  other  gas  sources  is  within  the  experimental 
scatter.  It  is  especially  easy  to  vary  pressure  in  the  shock  tube  at  high 
temperatures.  The  major  drawbacks  are  the  impulse  the  ultrasonic 
transducers  receive  from  the  shock  front,  and  the  presence  of  acoustic 
noise  from  boundary  layer  interaction  (Ref.  109)  and  discontinuities  in 
the  shock  tube.  At  present  iiltrasonic  measurements  in  the  shock  tube 
are  restricted  to  test  gases  with  at  least  50%  inert  gas  because  of  the 
noise  associated  with  boundary  layer  interaction  and  separation  in  poly¬ 
atomic  gases. 

* 
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3.6  ULTRASONIC  MEASUREMENTS  OF  THERMODYNAMIC 
PROPERTIES 

The  ultrasonic  velocity  may  be  used  to  obtain  the  equation  of  state 
of  moderately  dense  gases.  The  velocity  of  sound  measured  as  a  function 
of  temperature  and  pressure  is  sufficient  to  determine  the  compressibility 
if  independent  values  of  C  and  C  are  known  (Ref.  110).  The  sound  speed 
is  given  by  P  . 


Substituting  the  usual  equation  of  state 

p  =  zp  RT  (52) 

into' (51)  and  integrating,  we  obtain 
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The  experimental  values  of  7  (from  Cp  and  C^,)  and  c  as  a  function  of 
pressure  may  be  numerically  integrated  to  obtain  z  from  Eq.'  (53). 


When  experimental  determinations  of  y  are  not  available  from 
other  sources  a  theoretical  value  of  y  may  be  used  as  a  first  guess.  The 
estimated  7  (p)  is  then  used  to  obtain  a  first  approximation  to  z(T)  using 
Eq.  (53).  The  resulting  z(T)  then  may  be  used  to  obtain  improved  values 
of  7  from  the  thermodynamic  formulas  (Ref.  Ill)  for  Cp  and  ; 


C  -C 

V  V 


f 

•'0 


2pR  +  pRT  ^ 
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o  o 
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p  p  V  V 
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o  o 

where  and  are  the  ideal  gas  specific  heats.  This  process  m.ay 
be  repeated  until  the  best  approximations  to  7  and  z  are  obtained. 

Sound  speed  has  been  used  for  thermodynamic  property  measure¬ 
ments  by  many  investigators  (Ref.  55).  Using  the  technique  available 
at  Panametrics  the  thermodynamic  properties  of  gases  may  be  measured 
up  to  about  2000  K  at  pressures  up  to  about  10^  atm. 

A  typical  experimental  apparatus  which  has  been  used  up  to  about 
100  atm  is  shown  in  Figure  32.  Some  room  temperature  measurements 
taken  in  this  apparatus  are  shown  in  Figure  33.  The  difference  between 
the  increasing  and  decreasing  pressure  data  points  is  due  to  systematic 
deviation  in  temperature  due  to  compression  and  expansion.  The  tem¬ 
perature  difference  is  about  0.  8  K,  The  sound  speed  measurements  are 
compared  to  recent  measurements  by  other  investigators  in  Figure  34. 
The  sound  speed  measured  by  El- Hakeem  (Ref.  112)  is  lower  than  the 
Panametrics  measurements  by  an  amount  which  is  roughly  what  one 
would  expect  based  on  the  temperature  difference  of  the  measurements 
and  the  rule  of  thumb  c  ~  n/  7RT  • 


3.7  LASER  ACOUSTICS 


When  electromagnetic  radiation  passes  through  a  medium  some  of 
the  radiation  is  scattered  as  a  result  of  fluctuations  in  local  index  of  re¬ 
fraction  or  more  properly  dielectric  constant.  The  fluctuations  (Ref. 
113)  in  dielectric  constant,  £  ,  are  responsible  for  the  scattering  of  ra¬ 
diation.  These  fluctuations  are  caused,  in  turn,  by  the  fluctuations  in 
density,  temperature,  and  orientation  of  the  molecules.  In  gases  near 
one  atmosphere  the  fluctuations  in  the  number  density,  p  ,  completely 
predominate,  and  the  cross  section,  a,  per  unit  solid  angle  and  unit 
volume  for  polarized  light  can  be  shown  to  be  (Refs.  114,  115) 


(T  (K,  S2) 


X 


d^r 


dt  Ap  (?,  t)  Ap  (0,  0) 


X  exp  i  (K  •  r  -  fit)  , 


(56)' 
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where  O/Ztt  is  the  frequency  shift  of  the  scattered  radiation,  0  is  the 
angle  between  the  electric  field  of  the  incident  radiation  andtthe  wave 
vector  Kg  of  the  scattered  radiation,  is  the  wave  vector ]bf  the  in¬ 
cident  radiation  in  the  medium,  and  Ap  is  the  fluctuation  of/the  number 
density  about  its  average  p  q.  The  wave  vector  of  the  fluctuation  being 
observed,  K,  is  the  vector  difference  between  Kg  and  and  its  mag¬ 
nitude  is  given  by  K  =  2Kq  sin  (0  /2),  since  Kg  »  Kq  .  The  correlation 
function  for  the  fluctuations  in  density  is  related  to  the  classical  limit 
of  Van  Hove' s  density-density  correlation  function  (Ref.  Il6),  G(r,  t), 
by 


<  Ap  (?;  t)  Ap  (0,  0)  >  =  p^  (G(r,  z)  -  p  (57) 


Van  Leeuwen  and  Yip  (Ref.  117)  have  shown  that  for  a  dilute  gas 

<  Ap  (?,  t)  Ap  (0,  0)  >  =  p  ^  J  d%  f  (?,  p,  t)  (58) 

where  f  (?,  ^,  t)  is  the  particle  distribution  function  satisfying  the  lin¬ 
earized  Boltzmann  equation  subject  to  the  initial  condition 

f  (r;p,o)  =  (2Tr  mkT)“^/^  exp  mk  T)  6  (r).  (59) 

Here  mis  the  mass  of  the  molecule,  k  is  the  Boltzmann  constant,  and  p 
is  the  particle  momentum.  It  follows  that  O’  (K,  f2)  is  proportional  to  the 
double  Fourier  transform  of  the  space  and  time  distribution  function  for 
the  gas  molecules.  In  the  laser  scattering  experiment  the  scattered 
intensity  vs  wavelength  is  interferometrically  determined,  where  the 
trace  that  is  recorded  is  the  convolution  of  O’  (K,  O)  with  the  instrument 
profile  of  the  interferometer. 

The  form  of  the  spectrum  of  the  scattered  light  is  determined  by  a 
parameter  (Ref,  118)  y,  which  is  a  measure  of  the  ratio  of  the  wave¬ 
length  of  the  fluctuation  observed  by  the  collision  mean  free  path 

y  =  (2ir/K)  /  (27r/a)  (2kT/m)^/^,  (60) 
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where  a  is  an  effective  collision  frequency.  When  y  >  10  the  hydrody¬ 
namic  description  is  very  good.  The  spectrum  consists  of  three  distinct, 
nonoverlapping  lines;  an  unshifted  component  due  to  nonpropagating 
fluctuations  (the  Rayleigh  component)  and  a  symmetrically  shifted  doublet 
due  to  the  scattering  from  thermally  excited  sound  waves  (Brillouin 
scattering).  The  frequency  of  co  of  the  Doppler  shifted  Brillouin  com¬ 
ponents  is  related  to  the  sound  velocity  by 


c  =  CO  /K  . 


(61) 


As  a  result  of  acoustic  attenuation  the  Brillouin  lines  are  broadened 
and  assume  a  Lorentzian  shape  with  half  intensity  widths 


AcOb 


2  TK^  . 


(62) 


T  is  proportional  to  the  sound  attenuation  coefficient  and  is  given  as 


-  1  (3  n  +  T7*  .  \ 

-  2  \  p  pCp  -pC^/- 


(63) 


Since  the  time  decay  of  thermal  fluctuations  at  constant  pressure  is  de¬ 
termined  by  the  thermal  diffusivity  \  /p  Cp,  the  Rayleigh  line  is  also 
broadened  and  assumes  a  Lorentzian  shape  with  half  intensity  width 
Aco  given  by 

£x 


Aco, 


2\ 

pc 


K 


(64) 


Theoretical  considerations  from  the  thermodynamics  of  fluctuations 
(Ref.  113)  relate  the  area  under  the  Rayleigh  component  to  the  total 
area  under  the  Brillouin  lines  to  the  specific  heats  ratio,  C  /C^.  How¬ 
ever,  since  the  central  line  and  the  Brillouin  lines  all  have  Lorentzian 
shape,  the  area  ratio  becomes  just  the  ratio  of  intensities,  1^/2  1^. 
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(65) 
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the  Landau- Placzek  ratio  (Ref.  119)  for  intensities.  Figure  35  illus¬ 
trates  the  configuration  of  the  scattered  intensity  spectrum. 

In  the  region  2  <  y  <  10,  the  three  lines  begin  to  overlap  since 
their  widths  increase  faster  with  K  than  does  their  splitting;  yet,  the 
spectrum  can  still  be  described  rather  accurately  by  the  hydrodynamic 
equations.  However,  when  y  <  2  the  hydrodynamic  equations  no  longer 
adequately  describe  the  spectrum.  A  kinetic  theory  must  now  be  used  to 
take  into  account  the  effects  of  the  distribution  in  molecular  velocities. 

It  will  prove  to  be  instructive  to  examine  a  model  laser  gas  scat¬ 
tering  system  such  as  the  one  illustrated  in  Figure  36.  As  will  be  shown, 
the  wavelength  of  incident  light  must  be  strongly  monochromatic.  The  re¬ 
quired  high  spectral  purity  can  be  obtained  from  a  single  mode  stabilized 
laser.  Light  from  the  laser  scattered  in  the  scattering  cell  is  collimated 
by  a  conical  lens,  however,  since  the  scattering  amplitude  is  small  the 
unscattered  light  is  occluded  by  an  aperture  in  the  diagonal  mirror.  The 
unscattered  light  is  detected  and  fed  back  to  a  servo  control  system  for 
the  single  mode  laser.  The  scattered  light  is  reflected  by  the  diagonal 
mirror  into  a  sweeping  interferometer,  phototube  detector  and  pulse 
amplifier- discriminator. 

In  order  to  gain  some  notion  as  to  the  refinements  necessary  for 
the  experimental  apparatus  consider  the  test  gas  to  be  argon  at  NTP  and 
assume  for  purposes  of  simplicity  that  y  =  10.  Argon  at  NTP  has  a  mean 
free  path  of  10“^  cm  so  that  for  y  =  10  the  fluctuation  wavelength  becomes 

10“'^  cm.  Now  K  =  sin  —  ,  giving 

^  o  2 


sin  ^  =  .  3165  (66) 

for  the  He-Ne  laser  line  at  6330A  and 

9  =  36.9°  (67) 

as  the  .scattering  angle.  S2  =  2.  976x  10^^  radians/sec  corresponding  to 

\  =  633 OA.  From  co  =  Kc  for  the  Brillouin  line  shift 

o 

Q 

w  =  2.  03  X  10  radians/sec  .  (68) 
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The  ratio  of  center  frequency  to  the  frequency  shift 


£2 

CO 


1.  5  X  10 


6 


(69) 


In  order  to  unambiguously  detect  frequency  shift  of  these  magnitudes 
optical  instruments  with  resolving  powers  of  at  least  10  —  are  required, 

say,  of  the  order  of  2  x  10^.  This  high  degree  of  optical  resolution  is 
obtained,  for  example,  in  Fabry-Perot  interferometers. 

Backs catte ring  relaxes  somewhat  the  resolving  power  requirement 
since  for  backscatter  K  =  coK^  .  Here  y  =  3,  however,  so  that  analysis 
of  scattering  intensity  may  well  require  recourse  to  kinetic  theory  cal¬ 
culations.  cj  =  2  K  c  gives 


CO 


6.  41  X  10*^ 


(70) 


which  relaxes  the  resolving  power  requirement  by  a  factor  of  3, 

The  laser  acoustic  method  requires  sensitive  optical  and  electronic 
measurements.  Since  optics  must  be  kept  free  from  vibration,  elaborate 
precautions  may  be  required.  The  scattered  light  intensity  is  very  low 
and  thus  long  time  counting  is  required  to  extract  the  scattered  signal 
from  the  background.  As  the  gas  density'is  increased  the  amount  of 
scattered  light  increases.  From  this  point  of  view  the  experiment  be¬ 
comes  easier  to  perform  as  the  gas  density  is  raised  (see  Equations  56 
and  57). 

As  pointed  out  by  several  investigators,  the  laser  acoustic  tech¬ 
nique  may  be  used  to  measure  sound  speed,  ultrasonic  absorption  and 
thermal  diffusivity.  This  means  that  both  viscosity  and  thermal  con¬ 
ductivity  are  completely  determined.  With  the  usual  ultrasonic  tech¬ 
nique  only  the  sum  of  these  properties' is  determined.  In  addition,  if 
the  gas  temperature  is  known,  the  sound  speed  provides  a  check  on  the 
thermodynamic  properties  of  the  gas  (i.  e.  ,  the  compressibility).  In¬ 
ternal  modes  which  are  not  coupled  to  translation  (slow  modes  such  as 
vibration)  are  easily  frozen  out  because  with  visible  light  and  at  the 
smallest  practical  scattering  angles  the  ultrasonic  frequency  is  of  the 
order  of  100  MHz.  Since  laser  acoustic  technique  provides  information 
concerning  the  fluctuations  of  density  in  a  temperature  in  the  gas,  this 
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may  be  the  most  useful  result  to  be  obtained  from  the  laser  scattering 
experiment  (Ref.  120). 

3.8  ULTRASONIC  ABSORPTION  MEASUREMENTS  IN 
HELIUM 

The  acoustic  experiments  were  performed  during  the  past  year  in 
the  reflected  region  of  the  shock  tube.  The  temperature  range  of  the 
experiments  was  from  500  K  to  3000  K  at  pressures,  of  l/2  to  2  atmo¬ 
spheres.  Frequencies  of  1  to  2  MHz  were  employed.  The  measure¬ 
ments  were  made  both  in  the  2-l/2"  by  1-3/4"  shock  tube  and  a  boundary 
layer  stripper  which  gave  1"  and  1-1/2"  acoustic  path  lengths.  The  basic 
measurement  consists  of  measuring  the  amplitude  across  two  different 
path  lengths  of  test  gas.  Assuming  the  boundary  layer  thickness  is  the 
same  on  each  of  the  side  walls  of  the  test  section,  a  two  path  measure¬ 
ment  allows  the  boundary  layer  to  be  canceled  out. 

A  diagram  of  the  ultrasonic  apparatus  is  shown  in  Figure  37.  A 
shock  wave  is  generated  in  10  to  20  mm  Hg  of  helium  using  400  to  500 
psi  of  helium  or  hydrogen  as  the  driver  gas.  The  time  of  arrival  of  the 
shock  front  is  measured  by  two  side  wall  heat  transfer  gauges  and  an 
end  wall  heat  transfer  gauge.  The  first  side  wall  heat  transfer  gauge 
provides  a  reference  time  for  the  other  events  which  occur  later.  The 
signal  from  this  gauge  is  fed  into  a  three  channel  delay  generator.  One 
channel  triggers  the  dual  beam  oscilloscope  displaying  the  responses  of 
the  downstream  heat  transfer  gauges.  The  other  channels  trigger  the 
rf  oscillator  and  the  readout  oscilloscope  for  the  ultrasonics.  The  use 
of  delays  allows  an  expanded  time  base  to  be  used  on  the  oscilloscopes 
improving  the  accuracy  of  the  time  measurements.  The  trigger  of  the 
ultrasonic  oscillator  is  also  delayed  to  send  the  ultrasonic  pulse  through 
the  most  uniform  part  of  the  reflected  shock  region. 

There  are  several  calibrations  involved  in  the  data  reduction.  .  The 
time  bases  of  the  oscilloscope  are  calibrated  using  subharmonics  of  a 
local  radio  station  (WBZ).  The  amplitude  calibrations  (voltage  scale) 
are  accomplished  by  transmitting  pulses  through  a  known  gas,  such  as 
He  or  Ar,  at  room  temperature.  The  gas  and  the  pressure  are  chosen 
to  reproduce  the  signal  amplitude  obtained  in  the  high  temperature  gas 
as  closely  as  possible.  The  transmitted  amplitudes  in  the  300  K  tests 
calibrate  the  activity  of  the  ultrasonic  transducers.  The  sound  speed 
at  300  K  is  known  to  within  0.  1%.  Therefore,  the  transit  time  and  the 
measured  path  length  are  used  to  calibrate  out  any  systematic  delays  in 
the  transmit- receive  system. 
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The  frequencies  of  the  ultrasonic  wave  are  1  MHz  and  1.  89  MHz. 

The  sound  absorption  increases  with  temperature  and  follows  the  ab¬ 
sorption  calculated  from  the  theoretical  transport  data  of  Emmons  and 
Lick  rather  closely  (within  4%).  Since  the  Emmons  and  Lick  thermal 
conductivity  and  viscosity  are  based  on  the  most  recent  He-He  inter¬ 
action  potentials,  these  measurements  confirm  the  high  cross  sections 
at  close  interaction  distances. 

The  reduction  of  the  ultrasonic  measurements  to  individual  trans¬ 
port  properties  may  be  accomplished  by  several  methods.  The  optimum 
technique  depends  on  the  type  of  auxiliary  information  that  is  available. 

In  the  present  case  a  Eucken  factor  was  used  based  on  the  experimental 
data  and  theory.  At  the  lower  temperature,  T  <  2000  C,  the  viscosity 
measured  by  other  techniques  could  have  been  used  to  obtain  thermal 
conductivities  to  clear  up  the  large  discrepancy  between  various  mea¬ 
sured  values.  However,  since  the  contribution  of  thermal  conductivity 
to  the  absorption  is  only  30%  of  the  contribution  of  viscosity,  the  absorp¬ 
tion  data  were  reduced  to  viscosity.  The  viscosity  vs  temperature  is 
shown  in  Figure  38. 

The  viscosity  from  several  sources  is  shown  in  Figure  34  in  addi¬ 
tion  to  the  ultrasonic  measurements.  Both  of  the  measured  values 
(Guevara  and  Wageman,  authors)  fall  close  to  the  theoretical  values  of 
Lick  and  Emmons  indicated  by  the  solid  line.  Earlier  ultrasonic  mea- 
surements  up  to  1300  K,  which  were  obtained  in  a  muffle  tube,  are 
systematically  higher  than  the  present  shock  tube  measurements.  The 
calculated  viscosities  of  Amdur  and  Mason  are  within  a  few  percent  up 
to  1500  K.  However,  at  2500  K  the  Amdur  and  Mason  viscosities  are 
~16%  higher  than  the  more  recent  viscosity  values  indicated  on  the  graph. 
There  are  indications  in  Figure  34  that  the  Lick  and  Emmons  values  are 
too  low  by  about  2%  over  much  of  the  temperature  range  indicated  in 
Figure  24  and  that  the  discrepancy  increases  at  the  higher  temperatures. 

Since  helitun  is  a  monatomic  gas,  a  potential  which  gives  low  vis¬ 
cosity  values  also  predicts  low  thermal  conductivity  values.  The  ultra¬ 
sonic  technique  is  very  sensitive  to  this  type  of  error.  The  absorption 
due  to  thermal  conduction  must  be  subtracted  from  the  total  absorption 
to  obtain  the  viscous  losses.  Therefore  the  disagreement  between 
theory  and  experiment  is  approximately  doubled.  The  comparison 
between  theory  and  experiment  is  more  obvious  in  the  deviation  charge 
for  the  viscosity  of  helium  (Figure  1).  The  measured  viscosities  are 
4%  higher  than  theory  at  200  K,  within  1%  at  intermediate  temperatures 
~500  K,  and  appear  to  systematically  deviate  at  the  higher  temperatures. 
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Similar  deviations  are  apparent  (Figure  2)  between  the  theory  of  Lick 
and  Emmons  and  the  measured  values  of  thermal  conductivity  above 
lOOO^K.  Extrapolation  of  the  deviations  indicated  on  the  charts  to 
6000  K  would  put  the  transport  properties  of  Lick  and  Emmons  about 
15%  lower  than  experiment.  Further  measurements  using  the  heat 
transfer  and  ultrasonic  techniques  could  test  the  validity  of  this  small 
uncertainty.  At  any  rate,  the  above  discussion  illustrates  how  errors 
can  accumulate  when  potential  functions  are  used  to  extrapolate  trans¬ 
port  properties  much  beyond  the  range  of  experimental  verification. 
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APPENDIX  I 

MATHEMATICAL  RELATIONS  FOR  THE  END  WALL 
BOUNDARY  LAYER 

The  nature  of  the  end  wall  boundary  layer  may  be  seen  by  a 
closer  examination  of  the  energy  equation. 


8  /pX.  8h  \ 

dn  |Cp  hr)j 


.  8h 


=  0 


where 


rj  = 


^  f 


p  dy 


In 


1(4  (41=  ‘[(s).  (4l-/‘ 


(71) 


(72) 


o  \  0  h 

Dividing  each  term  by  ^ —  — —  and  integrating  between  the  wall  rj  =  0, 

and  some  point  at  the  outer  edge  of  the  boundary  layer  rj  j  we  get 


^  T)  dn.  (73) 


Since 


9n 


0  as  17  — >  00,  the  first  term  on  the  right  is  a  large  neg¬ 


ative  constant  if  we  chose  rj  j  to  be  a  point  near  the  free  stream.  The 
second  term  on  the  right  is  positive  and  depends  on  all. of  the  values  of 
for  all  the  temperatures  in  the  boundary  layer.  j  varies  al- 

irSst  linearly  with  rj  until  the  edge  of  the^boundary  layer  is  reached. 

At  this  point  the  goes  to  zero  as  e“*  .  Once  becomes  very 

small  the  integral  on  the  right  hand  of  (73)  sin^ply  cancels  the  first  term 
without  adding  significant  information  about  —  .  Therefore,  the  en¬ 
thalpy  gradient  at  the  wall  contains  significant^information  concerning 
the  thermal  conductivity  for  temperature  somewhat  less  than  the  free 
streeim  temperature. 


Equation  (73)  may  be  used  to  obtain  a  formula  for  the  thermal 
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conductivity  in  the  outer  edge  of  the  boundary  layer  in  terms  of  the  den¬ 
sity  gradients.  This  formula  is  useful  for  obtaining  thermal  conductivity 
from  optical  measurements  on  the  end  wall  of  the  shock  tube.  A  second 
equation  such  as  (73)  may  be  obtained  by  integrating  to  a  second  point 
near  the  outer  edge  of  the  boundary  layer.  Subtracting  these  equations, 
one  obtains 


In 


(74) 


Sufficiently  near  the  outer  edge  the  density,  specific  heat  and  thermal 
conductivity,  are  essentially  the  free  stream  values.  In  addition,  if 
density  is  known  as  a  function  of  enthalpy  (the  pressure  in  the  boundary 
layer  is  constant)  h  =  h  (p  )  then 


8h  _  a/TT  9h  _  sTzT  /8h  \  ££ 
drj  p  9y”  p  \9p/_9y 

e  p 


(75) 


The  new  thermodynamic  variable 

free  stream  density.  Under  these 
hand  side  of  (73)  becomes 


j  may  also  be  evaluated  at  the 
conditions  the  integral  on  the  right 


J, 

2 


C 

_2£ 

p\ 


and  may  he  expressed 


=  + 


V2t 


(72 -Yi)  • 


(76) 


(77) 


Finally,  Eq.  (74)  becomes 

■  (ty)  "  *‘’'2  ■  ’'l’^  ^  ®’’'2  ■ 


(78) 
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w.here  ■ 
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1 
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C 

X 

e 


8t 


(79) 


and  B  and  C  are  constants  which  follow  by  substitution  of  (76)  and  (77)  in 
(74).  Thus,  fitting  the  logarithm  of  the  measured  density  gradient  (from 
a  'schlieren  photograph)  to  a  quadratic  function  of  distance  such  as  (78), 
the  thermal  conductivity  may  be  obtained  from 


X 

e 


1 

2 


(80) 


A  different  approach  to  Eq.  (71)  leads  to  a  formula  which  may 
be  used  to  calculate  thermal  conductivity  from  heat  transfer  measure¬ 
ments.  The  independent  variable  in  (71)  may  be  changed  to  a  heat  flux 
potential 


dh  or  6*  = 

C  ^  C 

P  P 


drj 


(81) 


The  energy  equation  becomes 

C 

<(>•  +  n  ^  =  0.  (82) 


Integrating  (82)  between  zero  and  infinite  and  noting  that  =  0  at  in¬ 
finity,  one  obtains 


Using  (81),  this  becomes 


w 


r)  dh. 


(83) 


(84) 
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The  variable  r}  may  be  expressed  as  a  power  series  in  (p  using  a  Taylor 
expansion  about  ^  =  0.  The  coefficients  of  the  expansion  maybe  obtained 

by  calculating  from  (71).  Keeping  the  first  term  in  the  expression 

d</>" 

gives 


r? 


(85) 


or  substituting 

(85)  into  (84), 

integrating  the  result,  and  substituting  in 

(81)  we  obtain 

M 

m'  - 

h 

h‘ 

X 

^  dh  dh»  (86 

\  Vw 

w 

w 

P 

where  a  is  a  constant  which  corrects  the  integral  for  the  errors  due  to 

neglecting  higher  powers  of — ^ —  in  (85). 

6  * 

^  w 


Finally,  noting  that 


dn 


(87) 


near  the  wall  the  expression  for  the  heat  flux  is 
Q  _  ^  w  /dh  \  _ 

~  7W  ^  c  Idx/  = 

*^w 

Comparing  (88)  with  exact  numerical  integrations  of  the  energy  equation 
Kemp  showed  that  choosing  a  =  1.  13  gave  the  correct  answer  to  within 
1%  for  the  cases  of  interest  here.  In  practice  the  approximation  (the 
value  of  a)  (88)  should  be  checked  by  numerical  integration  of  an  ideal 
gas  similar  to  the  gas  under  investigation. 
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■  The  first  fact  to  notice  about  (88)  is  that  an  error  of  8%  in  the  mea¬ 
surement  of  Q  reflects  as  an  error  of  16%  in  the  integral 

h^  h' 

I  ^  ^ 

w  w  P 

Suppose  the  thermal  conductivity  of  the  gas  is  known  over  half  the  range 
of  integration.  Then  only  about  50%  of  the  integral  I  bears  significant 
information  on  thermal  conductivity  in  the  temperature  range  of  the 
outer  half  of  the  boundary  layer.  Therefore  the  uncertainties  in  the 
experimental  determination  of  I  are  doubled  in  the  thermal  conductivity. 
Similarly,  large  uncertainties  in  the  thermal  conductivity  over  a  small 
temperature  range  will  not  be  detected  in  the  heat  transfer  rate  to  the 
end  wall  of  the  shock  tube.  However,  systematic  discrepancies  over 
large  temperature  ranges  such  as  the  disagreement  between  thermal 
conductivity  of  helium  from  different  sources  (discussed  in  the  data 
review  section)  will  show  up  in  the  heat  transfer  measurement. 


69 


DEVIATION  IN  PERCENT 


8 


200  500  1000  1500  2000  2500 


TEMPERATURE  "K 


Fig.  1  Viscosity  deviation  chart  for  helium 
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Fig.  2  Thermal  conductiviry  deviation  chart  for  helium 
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EXCESS  VISCOSITY  POISE 


Fig.  3  Excess  viscosity  for  helium 
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Fig.  5  Viscosity  deviation  chart  for  argon 
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Fig.  6  Thermal  conductivity  deviation  chart  for  argon 
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Fig.  7  Thermal  conductivity  for  orgon  in  the  range  3000° K-1 2,000° K 
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Fig.  8  Excess  viscosity  for  argon 
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Fig.  9  Excess  thermal 
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Fig.  11  Viscosity  deviation  chart  for  hydrogen 
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Fig.  12  Thermal  conductivity  deviation  chart  for  hydrogen 
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Fig.  13  Eucken  factor  for  hydrogen 


AEDC.TR-68-105 


00 

CO 


6H 


Ul  2 

O 

(T 

Ul 

O. 

?  0 

z 

o 

p 

< 

>^2 

Ul 

a 


-6 


wA  A 


O  0  o  o. 


A 

□  □ 


iP— a 


'u- 


a^L 


□ 


□  A 
Q 


A« 


X 


X 


X 


REFERENCES 

Q132 

0150 

A151 

0130 


X 


1.295 


2.199 


2.913 


3.490  4.00  4.454  6.31 

REFERENCE  VALUES  XilO^  (dyne'-SEC/CM^ ) 


± 


X 


200 


500 


1000  1500 

TEMPERATURE  "K 


2000 


2500 


Fig.  14  Viscosity  deviation  chart  for  nitrogen 
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Fig.  16  Thermal  conductivity  deviation  chart  for  nitrogen 
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Fig.  17  Eucken  factor  for  nitrogen 
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Fig.  18  High  temperature  thermal  conductivity  deviation  chart  for  nitrogen 
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Fig.  21  Thermal  conductivity  deviation  chart  for  air 
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Fig.  23  High  temperature  thermal  conductivity  deviation  chart  for  air 
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Fig.  24 


INSULATION 


schematic  of  a  high  temperature  capillary  viscometer 
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Fig.  27  The  effect,  of  ore  current  and  diameter  of  the  effective  thermal  conductivity 
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REFLECTED  SHOCK  REGION 
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Fig.  28  Heat  transfer  in  the  reflected  shock  region 
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Fig.  29  Experimental  apparatus  for  schlieren  measurements 
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Fig.  30  Schematic  of  schlieren  system 


.  31  Schlieren  photograph  of  end  wall  boundary  layer 
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Fig.  32  Test  chamber  for  high  pressure  sound  speed  measurements 
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Fig.  34  Velocity  of  sound  in  nitrogon  gas  os  a  functionof  pressure 
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Fig.  35  Scattered  intensity  spectrum 
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Fig.  36  Model  loser  gas  scattering  system 
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Fig.  37  Diagram  of  apparatus  for  ultrasonic  measurements  in  the  shock  tube 
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Fig.  38  Experimental  viscosities  for  helium 
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